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ABSTRACT 

The magnetism and rotation of white dwarf (WD) stars are investigated in relation to a hydromag- 
netic dynamo operating in the progenitor during shell burning phases. The downward pumping of 
angular momentum in the convective envelope, in combination with the absorption of a planet or tidal 
spin-up from a binary companion, can trigger strong dynamo action near the core-envelope boundary. 
Several arguments point to the outer core as the source for a magnetic field in the WD remnant: the 
outer third of a ~ 0.55 Mq WD is processed during the shell burning phase(s) of the progenitor; the 
escape of magnetic hclicity through the envelope mediates the growth of (compensating) helicity in 
the core, as is needed to maintain a stable magnetic field in the remnant; and the intense radiation 
flux at the core boundary facilitates magnetic buoyancy within a relatively thick tachocline layer. 
The helicity flux into the growing core is driven by a dynamical imbalance with a latitude-dependent 
rotational stress. The magnetic field deposited in an isolated massive WD is concentrated in an outer 
shell of mass < 0.1 Mq and can reach ~ 10 MG. A buried toroidal field experiences moderate ohmic 
decay above an age ~ 0.3 Gyr, which may lead to growth or decay of the external magnetic field. The 
final WD spin period is related to a critical spin rate below which magnetic activity shuts off, and 
core and envelope decouple; it generally sits in the range of hours to days. WD periods ranging up to 
a year are possible if the envelope re-expands following a late thermal pulse. 

Subject headings: stars: giants - stars: white dwarfs - stars: rotation - magnetic fields 


1. INTRODUCTION 

The magnetic field and spin of a white dwarf (WD) 
are inherited from its progenitor star. Our interest here 
is in how magnetism and rotation evolve together within 
the progenitor, in response to an active hydromagnetic 
dynamo. We focus on the post-main sequence phase 
of stellar evolution, when the material comprising the 
outer layers of the WD is deposited by a strong inflow 
from a hydrogen-rich convective envelope, through burn¬ 
ing shell(s), into the stellar core. The presence or absence 
of a magnetic held at the surface of the WD remnant de¬ 
pends on the ability of the inner envelope to sustain a 
magnetic held. 

Rotationally-driven dynamo activity on the red gi¬ 
ant (RGB) and asymptotic giant (AGB) branch has re¬ 
ceived only limited at tention (e.g. iBlackman et al.ll2001b 
iNordhaus et all [200811 . due to a prevailing assumption 
that the radial angular velocity prohle fl(r) in the con¬ 
vective envelope will mirror that observed in the Sun. 
Although localized magnetic held amplification in the 
outer, slowly rotating envelope m a y result in weak sur - 
face magnetic helds dDorchl [2004 lAuriere et all 1201 011 . 
a dynamo that manifests itself through strong chromo¬ 
spheric and coronal emission appears to depend on an 
interactio n with a binary ste llar companion, as in RSCvn 
stars (e.g. [Moss et ahlll991lh Magnetic activity appears 
to be prevalent during the first stages of p ost-MS ex¬ 
pansi on, and during core helium burning (JAuriere et al.l 
120141 . but the sampling of more expanded giants is very 
limited. 

On the other hand, the angular velocity prohle in a 


giant is less sensitive to boundary conditions than it is 
in the Sun, given the extreme depth of the envelope. 
Here we are guided by i) recent measurements of rapid 
core rotation in subgiants and core He burni ng stars us¬ 
ing asteroseismic measur ements by Kepler (Beck et al.1 
l2012t IMosser et aI1l2012ll : ii) measurements of a strong 
decrease of surface rotat ion speeds in subgiant s as they 
evolve toward the RGB (jSchriiver fe PolsHl993T) ; and iii) 
anelastic simulations which show a tendency to nearly 
uniform specihc angular mo mentum in a deep conve ctive 
envelope with slow rotation (iBrun fe Palacios 2009 i. 

In a companion paper (jKissin fe Thompsonll2015l here- 
after paper I), we show that the Kepler measurements 
are consistent with a tight coupling between the radia¬ 
tive core and convective envelope, combined with some 
inward pumping of angular momentum deep into the en¬ 
velope. The favored rotation prohle has H(r) transition¬ 
ing between oc r _1 in the inner envelope to oc r~ 2 in the 
outer and more slowly rotating layers. The depth of the 
transition point depends on the size of the star and its 
angular momentum. 

An important consequence of such an inward-peaked 
rotation prohle is that the threshold for a rotationally 
driven dynamo may be reached in the inner envelope, 
but not near the stellar photosphere. A key hgure of 
merit is the Coriolis parameter 

Co — ^T'con — (1) 

^con 

as measured at the base of the envelope; here v con is 
the convective velocity and T con is the timescale for con- 
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vection over one pressure scale height. Observations of 
MS stars indicate the presence of dynamo activity when 
Co > 1 (c orr espon ding to a Rossby number < 27 t; e.g. 
iReiners et al.112009 b The radial transition in the slope 
of the rotation profile in a deep convective envelope is 
postulated to occur where Co approaches unity, with the 
most rapidly rotating zone lying below this transition. 

Inward pumping of angular momentum can also reduce 
the need for magnetically driven spindown in subgiants 
(Paper I). 

1.1. Remnant Magnetism: Dependence on Stellar Mass, 
Binarity and Age 

Three hypotheses have dominated previous discussions 
of the origin of magnetic fields in WD stars: the conserva¬ 
tion of mag netic flu x from the main seq uence (MS) pro¬ 
genit or (lAneel et all 119811 : iFerrario fe W ickramasinghc 
l2005ll : the amplification of such a seed field by con¬ 
vective episo des in the core of the progenitor star 
(iRuderman fc Sutherland! 1 973 ): and the transient am¬ 
plification of the field durin g a merger event (|Tout et al.l 
120081: iNordhaus et alJl2011l) . 

We begin our investigation by reviewing some con¬ 
straints on each of these hypotheses; and then summa¬ 
rize the arguments pointing to an amplification of the 
magnetic field during the post-MS expansion of the WD 
progenitor. This process comprises three steps: i) the 
advection of angular momentum inward through the con¬ 
vective envelope; ii) the amplification of a magnetic field 
by a dynamo near the core-envelope boundary; followed 
by iii) a flow of magnetic helicity into the growing core, 
which is sustained by a persistent angular velocity gra¬ 
dient and magnetic twist below the boundary. 

Whether the inner envelope can support a dynamo 
turns out to depend on a combination of stellar prop¬ 
erties. The core rotation rate is sensitive to the loss of 
angular momentum on the MS, and to any interaction 
with a planetary or stellar companion. Clearly solid ro¬ 
tation would imply Co <C 1 throughout the envelope of a 
star approaching the tip of the RGB or AGB; and there¬ 
fore suggest the absence of a large-scale hydromagnetic 
dynamo. But when the envelope rotation profile is in¬ 
wardly peaked, we find differing answers for isolated stars 
of initial mass AJzams less or greater than ~ 1.3 Mq, cor¬ 
responding to the presence or absence of a strong mag¬ 
netic wind on the MS. Hence the greater importance of 
binary interaction for sustaining Co > 1 in solar-mass 
stars. 

Another consideration is the proportion of the ma¬ 
terial contained in the remnant WD that is deposited 
during single and double shell burning. This comprises 
about a third of the mass of the WD when Mzams ~ 1 
(M w d ~ 0.55M 0 ), but is strongly reduced by dredge- 
up (after core helium ignitition and before the onset of 
thermal pulses) in the intermediate-mass stars that are 
the progenitors of massive WDs. In this second case, the 
magnetized mass shell turns out to be thin enough that 
ohmic decay is initiated on a timescale ~ 0.3 Gyr (with 
some dependence on the mass loss rate on the AGB). 

The source of a strong magnetic field (1 — 10 3 MG) 
in a WD therefore depends on both its age and mass. 
Indeed, expanded sampling of WD magnetism calls into 
question the flux-freezing hypothesis. Whereas between 
~ 2 and ~ 10 percent of isolated WDs are magnetic in 


this sense dLiebert et al.l 1200311200511 - roughly similar 
to the proportion of chemically peculiar and magnetic A 
stars - t h e proport ion is higher in isolated massive WDs 
(iKenler et, al.ll2013h and extends to ~ 2 5 percent in cata¬ 
clysmic variables (CVs: lGhnsickell2005lf . A large propor¬ 
tion of magnetic WDs in short-period CVs also appear 
to have very strong magnetic fields, exceeding 10 MG, 
which are much rarer in the isolated WD population. 

It has also been noted that magnetic WDs appear to 
be rare in binar y system s with a non-interacting MS 
companion (|Liebert et alJl2005h . suggesting that binary 
interaction may be involved in their genesis. Previ¬ 
ous theoretical ideas have centered around the impul- 
sive growth of a magnetic field during a binary merger 
(INordhaus et al.l [20071 : iTout et all 120081 INordhaus et al.1 
I2011H . Here our focus is on a much longer-lived process 
that is sustained by the inward pumping of angular mo¬ 
mentum. 

A solar-mass star is slowly rotating at the end of the 
MS, so that the absorption of a giant planet signifi¬ 
cantly spins up the star and may trigger dynamo activity 
dLivio fc Sokeil[2002} ). The critical planetary mass that 
has this effect near the base of the envelope depends on 
the angular velocity profile. For a flat distribution of spe¬ 
cific angular momentum at Co < 1, we find that the ab¬ 
sorption of even a Neptune-mass planet can significantly 
spin up the core. 

A star more massive than Mzams ~ 1.3 Mq retains 
enough angular momentum at the end of the MS that 
only a planet more massive than Jupiter will significantly 
supplement its rotational angular momentum. Nonethe¬ 
less we find that such an interaction can lead to stronger 
magnetic fields in the WD remnant. 

Two properties of giants are especially relevant here. 
The inward flow of matter across the core-envelope 
boundary allows a net transport of magnetic helicity, 
which has been shown to be a necessary ingredient in 
t he long-term sta bilization of white dwarf magnetic fields 
dBraitliwaite fe Spruitl I2004H ; and has also been impli¬ 
cated in the establishm ent o f solid rotat ion in the radia¬ 
tive core of the Sun ( Gough fe Mclntvreli19981 ). We show 
that this helicity flux can be related to the magnitude 
of the large-scale Maxwell torque that enforces (nearly) 
solid rotation in the outer core. 

The outer core of an AGB star is also distinguished 
by an intense radiation flux. This greatly facilitates the 
buoyancy of a toroidal magnetic field that is anchored 
in a tachocline layer below the convective envelope. A 
thickness approaching ~ 10 percent of a scale height can 
maintain contact with the envelope over a dynamo pe¬ 
riod ~ 10 — 10 2 T con - This suggests that a given level of 
magnetic activity (measured e.g. by the ratio of Maxwell 
to Reynolds stresses) may be achieved with a lower Cori¬ 
olis parameter in an AGB envelope than in a solar-type 
MS star. Magnetic ally induced mixing o f core material 
into the envelope (INordhaus et al.l I2008D could also be 
significantly enhanced by this effect. 

1.2. Remnant Spin 

The measured spins of WDs are strongly inconsistent 
with solid rotation in the envelope of an AGB star, com¬ 
bined with a tight core-envelope coupling. But even if 
we allow for inward convective pumping of angular mo¬ 
mentum, the final spin rate of a WD remnant is sensitive 
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to the heavy loss of mass during the transition to the 
post-AGB phase]]] 

A decoupling of the rotation of core and envelope is 
only feasible when the large-scale poloidal magnetic field 
drops below the equivalent of 10 2 — 10 3 G in the WD 
during the brief, thermally pulsating AGB phase. (Much 
weaker fields will couple the core and envelope of a sub¬ 
giant or core helium burning star; Paper I.) This appears, 
at first sight, to be inconsistent with strong magnetism in 
the WD. In fact, a magnetic coupling between core and 
envelope is naturally self-limiting: once Co drops below 
a critical value, the dynamo operating in the inner enve¬ 
lope shuts off. A reasonable first guess (Co cr it ~ 0.1—0.3) 
yields WD rotation periods in the range observed in most 
WDs. 

As a result, the remnant rotation is inversely corre¬ 
lated with the net angular momentum in the star: the 
slower the rotation of the AGB envelope, the sooner the 
decoupling between core and envelope, when the core 
has a larger radius and moment of inertia. The weakly 
magnetized surface layer of the WD is very thin, and we 
show that the magnetic field below quickly diffuses to the 
photosphere. 

Our most surprising finding is that WDs with very 
long (~ yr) rotation periods originate from stars which 
spin relatively rapidly near the end of the AGB phase. 
Then the rotation of core and envelope remain magneti¬ 
cally coupled during the contraction of the envelope. The 
endpoint is sensitive to details of how the envelope is ex¬ 
pelled: in particular, to whether the envelope experiences 
a re-expansion following a late thermal pulse. This taps 
most of the remaining angular momentum of the core. 
Without such a re-expansion, the final rotation period 
can be shorter than a day. 

It should be kept in mind that anisotropic contraction 
of the e nvelope can contribute to white dwarf rotation 
ijSpruit] Il998f) . Indeed, it may be needed to generate 
centrifugally supported disks that provide seeds for metal 
enrichment of white dwarf atmospheres. 

1.3. Stellar Models and Planetary Interaction 

In order to evaluate various pieces of the physics, we 
employ two stellar models of solar metallicity and ini¬ 
tial mass Mzams = 1 M© and 5 AG Both mode ls are 
constructed using the MESA code dPaxton et alJ [20Tll . 
version 5527), with rotational degrees of freedom turned 
off. The normalization of the mass-loss rate on the AGB, 
and the masses of the WD remants, were fixed by match¬ 
ing the latest initial- to-fi nal mass r elation obtained from 
open clusters dDobbie et all l2009lf . The coefficient en¬ 
tering into the formula of iBlockerl (|1995[ 1 was taken to 
be r]B = 0.05. The corresponding masses of the WD 
remnants are 0.55 Mg and 0.87 Mq. 

The rotational evolution was followed in post¬ 
processing using the rotation models described in Sec¬ 
tion 12.11 The loss of angular momentum due to mass 
loss on the RGB and AGB was taken into account, along 
with the interaction with a companion planet of a range 
of masses: Earth, Neptune, Jupiter in the 1 M 0 model, 
and even heavier companions in the 5 Mq model. In 
the first case, the companions were started at a range 
of semi-major axes and orbital eccentricities. The orbit 

1 We thank Peter Goldreich for emphasizing this point to us. 


of the planet, and the exchange of angular momentum 
with the star, were computed taking into account tidal 
drag and the gravitational quadrupole that is generated 
by convection in the giant envelope. Further details are 
given in Paper I. 

1.4. Plan of the Paper 

Section [5] and [3] discuss several processes that are rele¬ 
vant to the growth and maintenance of a stable magnetic 
field in the core of a post-MS star, and which determine 
the spin angular momentum trapped there at the end of 
the AGB. Section [4] presents detailed estimates of WD 
dipole magnetic fields. Section [5] analyzes the magnetic 
decoupling of core and envelope near the tip of the AGB, 
and provides estimates of the final WD spin. Ohmic 
transport through a thin, unmagnetized surface layer in 
a WD, and the thicker magnetized layer below it, is con¬ 
sidered in Section [6] 

Various outstanding issues - such as the sensitivity of 
our results to the assumed angular velocity profile in the 
progenitor envelope, the impulsive growth of a magnetic 
field during a merger, and effects associated with the 
anisotropic contraction of the residual AGB envelope - 
are addressed in the concluding Section]?] The Appendix 
gives further details of the accumulation of magnetic he- 
licity in the stellar core. 

2. ROTATION AND MAGNETISM IN A RAPIDLY 
EVOLVING GIANT STAR 

Here we summarize the rotation model for giant stars 
that was developed in Paper I. We then review some 
physical constraints on generating sustainable white 
dwarf magnetic fields, with a focus on the accumula¬ 
tion of magnetic helicity during the rapid growth of the 
hydrogen-depleted core of a giant star. We also sum¬ 
marize the relative effectiveness of different convective 
episodes at sourcing WD magnetic fields. In this regard, 
there are some significant structural differences between 
the evolving cores of solar-mass and of intermediate-mass 
stars. 

2.1. Combined Effect of Convection and Poloidal 
Magnetic Field on Angular Momentum 
Redistribution 

Our understanding of the origin of WD rotation has 
been impeded by an incomplete knowledge of mixing pro¬ 
cesses in the radiative parts of the progenitor. Even a 
very weak poloidal magnetic field (corresponding to a 
magnetic field <C MG in the WD under the assumption 
of magnetic flux freezing) will enforce solid rotation in 
the evolving core. 

A popular approach recently has been to assume that 
such a large-scale poloidal field is initially absent, and to 
imagine that the radial magnetic field that is needed for 
angular momentum transport is generated by a current- 
drive n (‘Tay l er’) instability of a predom inantly toroidal 
field dSpruitl l2002t iCantiello et al.l 12014 and references 
therein). This leads to an incomplete coupling between 
core and envelope. We argue that even an inefficient dy¬ 
namo process operating near the core-envelope boundary 
will easily short-circuit this effect. 

We develop, as a working alternative, a simple and de¬ 
terministic model of angular momentum redistribution in 
evolving stars. This model applies to stars which retain 
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(or gain) enough angular momentum to support some 
dynamo activity near the core-envelope boundary. 

1. The outer core is assumed to rotate as a solid body, as 
enforced by large-scale Maxwell stresses, and to co-rotate 
on the average with the inner envelope. 

2. Inhomogeneous rotation is sustained mainly in con¬ 
vective regions of the star, especially the deep outer en¬ 
velope that forms following the completion of core H and 
He burning. 

The rotation profile that is formulated and calibrated 
in paper I can be summarized as follows. Where the 
Coriolis force can be neglected, outside some transition 
radius R c , the envelope maintains uniform specific angu¬ 
lar momentum, 

H(r) = Q(R C ) ; r > R c . (2) 

Here H is the angular velocity averaged over a spherical 
shell. We identify R c with the base of the envelope, 

Rc -^benv (^con ^ 1 Rhenv') i (H) 

when_the angular momentum of the star is small enough 
that H(r)r con < 1 throughout the envelope. Throughout 
this paper, the subscript ‘benv’ refers to the base of the 
envelope. 

When the star has more angular momentum, we main¬ 
tain the profile ([ 2 ]) in the outer part of the convection 
zone, and consider 


n(r) = fi(i2e) 

_ Cotrans / T_\ 
Tcon(Rc) \RcJ 


Rbenv ^ r <C R c (4) 


in the inner envelope. The Coriolis parameter at the 
transition radius Cotrans is taken to be unity, and the 
radial index 1 < a < 3/2. As the angular momentum 
increases R c moves outward, and may reach the surface of 
the star, in which case Co is larger than unity everywhere 
in the envelope. 

As we discuss in Paper I, a = 1 best matches the 
observed core rotation periods of sub-giant and helium 
burning stars, which we use to calibrate our rotational 
model. However, in stars of much greater size, such as 
those approaching the tips of the RGB and AGB, the 
core is more centrally concentrated and the gravitational 
field in the inner envelope, g(r) oc r~' 3 , steepens from 
/? = 1 to /3 = 2. Then the value of a that results from 
the inward transport of angular momentum by deep con¬ 
vective plumes increases from (l+/3)/2 = lto3/2. It 
is during this more advanced phase of evolution that the 
helicity flux into the core is expected to peak, and the 
rotational fate of the WD is determined. We therefore 
use a = 3/2 in our analysis. 


2 . 2 . Magnetic Helicity Growth in Giant Cores 

Net magnetic helicity H = J A ■ B dV appears to be 
an essential ingredient for maintaining long-term mag- 
netostatic stabi lity in a sta tic and entirely fluid star 
(lBraitliwaite fc Spruitl I2004T) . Here B is the magnetic 


field and A the vector potential from which it is derived. 
H is a topological charge in a perfectly conducting fluid 
without boundary. For example, when a closed poloidal 
magnetic loop carrying flux <f> p is also twisted through 
an angle ip, the helicity is H = = ^<f> p , where 

<f >0 is the toroidal flux. TL evolves only on a very long 
ohmic timescale in an extended astrophysical fluid un¬ 
less there is an exchange of magnetic twist with the ex¬ 
terior of the fluid, or with a medium of low electrical 
conductivity. 

The accumulation of helicity considered here is asso¬ 
ciated with rapid changes in the structure of the star. 
We expect this to be possible in the outer core and inner 
envelope of a giant star, since this zone is connected to 
the surface of the star by rapid convective motions. Core 
convection has a topology less favorable to the accumu¬ 
lation of H, as does the erasure of differential rotation 
by the stretching of magnetic field lines in the interior of 
a radiative star. A simplified time-dependent model of 
an i ntern al s hear-driven dyna mo has been developed by 
IWickramasinghe et al.1 (120141 . but does not address the 
growth of H. 

Such a secular accumulation of magnetic helicity - al¬ 
though occuring relatively rapidly near the tips of the 
RGB and AGB - is nonetheless a gradual process when 
measured over an individual dynamo cycle. There has 
also been longstanding interest in the role that helicity 
conservation may play in limiting the exponential growth 
of a seed magnetic field in a turbulent magnetofluid 
(IGruzinov fc Diamondl 11994 IBlackman fc Field! 120004 
[Vishniac fc Choi 120011 ). In mean-field dynamo models, 
this involves a suppression of the a effect, e.g. the 
toroidal electromotive force which sources a large-scale 
poloidal magnetic field. 

This limitation on fast dynamo action could be re¬ 
moved by the transport of magnetic helicity on a rela¬ 
tively short timescale, comparable to the rotation period. 
Such transport may in volve the removal of twist across 
the surface of the star (IBlackman fc FieTdl l2000bT) ■ or al¬ 
ternatively an internal flow of helicity betwe en rotational 
hemispheres (lYishniac fc Shapovalovll2014f) . In either of 
these cases, rapid dynamo growth is not directly depen¬ 
dent on changes in the net stellar helicity, obtained by 
summing the contribution from both hemispheres. The 
longer, secular effect being considered here necessarily 
depends on a breakdown of reflection symmetry between 
magnetic hemispheres; and, because it involves a change 
in an integral property of the star, must depend on the 
transport of twist across the stellar surface. 

The connection between helicity growth and mass in¬ 
flow can be examined with a simple model of a toroidal 
magnetic twist superposed on a radial magnetic field, all 
subjected to a radial MHD flow of speed v r . (The sign of 
this flow, relative to the convective boundary, fluctuates 
during thermal pulses on the AGB but is negative on the 
average.) 

Within a sphere coinciding with Rbenv, one finds 

Sv r J dS 

(5) 


d ^ = l[A.BdV 
dt dt J 


cm 

dt 


dRbe 


dt 


dSAm m = - 
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(Appendix®. The gauge 


base of the convective envelope, 


A<f,( 0 ,r) 


is convenient, where 


1 

r sin# 


$r(®) 


r e 

$ r (0) = / sin O'r 2 B r (9')d9' 

Jo 


( 6 ) 

(7) 


Bi 


47r p 



/Co r An 

\27 x ip n 


2 

con 


The Maxwell stress works out to 

BrBt/, prAr(An ) 2 

47t 2irN dyn 


( 11 ) 


( 12 ) 


is (27 t) _ 1 times the radial magnetic flux between the ro¬ 
tation axis and polar angle 9. Clearly $ r (7r) = 0. 

One notices that the surface integral appearing in 
Equation ([5]) has the dimensions of a torque: it is com¬ 
parable in magnitude to the Maxwell stress B r B<f,/4:ir , 
multiplied by the lever arm r sin 6, and integrated over 
the core-envelope interface. Even if the net torque acting 
on the core boundary vanishes, dB/dt may remain finite 
as a consequence of the different angular dependences of 
and B r r sin 9. 

We are interested here in the case where B r B^/47 t is 
the dominant stress driving the core material toward 
solid rotation. Then we can estimate the helicity flux 
into the core once we know the source of the Maxwell 
torque. The largest effect turns out to be the latitude 
dependence of the Reynolds stress that is applied to the 
core at the base of the convective envelope. 

We may estimate this stress by considering a dynamo 
process at work in a tachoclinc layer right below the 
convective boundary. Convection above the tachocline 
supports strong latitudinal variations in n which im¬ 
pose a radial mis-match Af2 with the mean core rotation 
f^benv = ff(-Rbenv)- The tachocline can therefore be iden¬ 
tified with the radial layer (of thickness A r) in which Af2 
is concentrated, 


vn ~ 



An, 

"Ar" ? 


( 8 ) 


Our focus here is on the parts of the tachocline where 
dn/dr > 0 (in the Sun, it is the equatorial band at 
latitudes < 45°). The magnetorotational instability 
(|Balbus fc Hawlevl I1994H is excited where dn/dr < 0, 
and is sustained by rapid radiative diffusion even in 
t he presence of sta ble radiative convective equilibrium 
(|Menou et al.ll2004l) . 

The poloidal magnetic field threading this part of the 
tachocline is wound up linearly on timescales short com¬ 
pared with the dynamo period P<jyn, which we normal¬ 
ize as AdynAfW 1 . We limit the magnetic field strength 
in the tachocline by noting that a radial magnetic field 
stronger than 

Br ~ (4 *p) 1/2 ^- (9) 

Tlyn 


will act to redistribute angular momentum in a layer of 
density p and thickness Ar, thereby shorting out the an¬ 
gular velocity gradient. The ratio of toroidal to poloidal 
held components is obtained from the linear winding 
term in the induction equation, 


27 t r An 

- B(k ~ - - - Br 

Pdyn * A r 


( 10 ) 


One finds that the toroidal magnetic field is in approx¬ 
imate equipartition with the rotational motions at the 


We emphasize that this estimate of the stress has been 
obtained from dynamical considerations, and side-steps 
the challenging question of how the a effect is actu¬ 
ally manifested in a dynamo operating near a radiative- 
convective boundary. For recent approaches to this prob¬ 
lem, in the context of lower-luminosity MS stars, see 
iBlackman fc Thomas! ( 20151 ) and references therein. 

We note that Afl ~ D is seen in anelastic calcula¬ 
tions of de ep and rapidly rotatin g (Co > 1) convective 
envelopes (jBrun fc Palacios! I2009H . and is suggested by 
the analytic treatment of angular momentum pumping 
in Paper I. The latitudinal gradient is milder in the Sun, 
where the equator rotates more rapidly than the poles, 
(n eq — flpoie )/ n ~ 0.1; but the Solar convective envelope 
also has a very different aspect ratio. 

Now consider the helicity flux integrated over one rota¬ 
tional hemisphere. Even in a state of rotational equilib¬ 
rium - where the angle-integrated Maxwell torque van¬ 
ishes - the integral (0 generally remains finite. Notice 
also that it is the product of B r and B^ that is regulated 
according to Equation m- The Maxwell stress does not 
depend on the sign of B r that is imposed by a dynamo 
operating near the core-envelope boundary, nor does it 
depend directly on the coherence of the magnetic field. 

The sum of the two hemispheric integrals for dB/dt 
would vanish if the magnetic field were reflection- 
symmetric about the rotational equator. But in gen¬ 
eral it is not: the distribution of radial magnetic flux 
across the northern hemisphere need not precisely mir¬ 
ror that threading the southern hemisphere. (Only the 
net spherical flux vanishes.) Such an incomplete hemi¬ 
spheric cancellation allows a net flow of magnetic twist 
into the core. 

Other sources of helicity flow can be considered. The 
inflow through the core-envelope boundary itself drives 
dn/dr, but on a much longer timescale than the rotation 
period. 

The WD magnetic fields that result from these two 
sources of B are compared in Section [3] 


2.3. Relative Importance of 
Different Convective Episodes 

The magnetic field of a WD is influenced by a se¬ 
ries of convective episodes in the progenitor. One mea¬ 
sure of the relative importance of different convective 
episodes is provided by the convective Mach number 
■Adcon = v con /c s (here c s is the adiabatic speed of sound), 
which in a stellar core can be estimated as 


M 


con 


f RcoieVcon 
\ GMcore 


1/2 


(13) 


As the core continues to contract during later stages of 
nuclear burning, its magnetic energy and gravitational 
energy scale in the same way with central density. To 
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the extent that the generated Maxwell stress is limited by 
the convective stress, e.g. B$ < (47rp) 1//2 u con , episodes 
of greater A4 CO n are capable of generating stronger mag¬ 
netic fields. 

One encounters much larger Af C on on the giant 
branches (~ 10 -3 — 10” 2 ) than during steady core H or 
He burning in intermediate-mass stars (~ 10 -5 — 10 -4 ). 
It is also worth noting that the intense period of convec¬ 
tion during the AGB phase is also the shortest in dura¬ 
tion. That implies a weaker cancellation in the magnetic 
helicity that accumulates in the core, due to a shifting 
distribution of poloidal magnetic flux (see Section 12.21) . 

Magnetic fields that are generated during a certain 
convective phase will remain trapped unless exposed to 
later stages of convection, or the exterior mass shells are 
expelled. As we examine in Section [6l ohmic diffusion 
in the WD remnants of intermediate-mass stars is re¬ 
stricted to a relatively thin magnetized layer. The rele¬ 
vance of dynamo activity on the RGB, horizontal branch 
and AGB is determined, in good part, by the maximum 
inward penetration of the convective envelope: one re¬ 
quires M(i?benv) A Mwd- Only in relatively massive pro¬ 
genitors (Mzams A) 4 M 0 ) does MS core convection ex¬ 
tend to the mass boundary of the WD remnant, where 
it can influence the visible magnetic field; we ignore this 
effect here. 

The core helium flash occuring in stars of mass < 
2.3 Mq is confined to a central subset of the eventual 
WD material. Because a star does not bec ome fully con¬ 
vective during the flash (e.g. iBildsten et alJl20l2il . the 
magnetic helicity of the convecting core material is essen¬ 
tially conserved. A similar conclusion applies to the brief 
convective layers that are triggered by thermal pulses in 
the helium-burning shell near the tip of the AGB. 

2.3.1. Stars with Minimal Dredge-up of Helium 
(Mzams < 2.3 Mq) 

In stars of initial mass Mzams ~ 1 Mq, magnetic fields 
generated on the RGB are buried inside the WD remnant 
under a layer of CO-rich material that is produced dur¬ 
ing double shell burning (Figure Q]). They are ohmically 
coupled on a timescale ~ 10 s — 10 9 yr to a thinner mag¬ 
netized shell that appears on the AGB (Section O. 

The ohmic timescale across the entire magnetized shell 
is longer than ~ 10 9 yr, meaning that any magnetic field 
deposited near the center of the star during the pre-MS 
phase would influence the surface of the star only on a 
timescale of a few Gyr. 

We conclude that the existence of young magnetized 
WDs with masses < 0.6 Mq points to dynamo activity 
in the convective envelope on the RGB and AGB. 

2.3.2. Intermediate-mass Stars (Mzams > 2.3 Mq ) 

Stars of intermediate mass are distinguished from those 
of solar mass in that a dynamo operating on the RGB 
cannot contribute directly to the WD magnetic field: the 
outer parts of the helium core are dispersed by dredge 
up before the C/O core is assembled. In addition, only 
a limited mass is converted to carbon and oxygen by 
double shell burning (Figure [2j). The fraction of the WD 
mass contained by this outer shell is sensitive to the rate 
of mass loss on the thermally pulsating AGB. 

The ohmic conduction time across this outer magne¬ 
tized shell is several 10 s years. The surface magnetic field 



1 - !io .9 (years) 



t - t 109 (years) 


Fig. 1 . — Interior profile of Mzams — IMq star around the tip 
of the RGB and first stages of core He burning (top panel) and dur¬ 
ing AGB phase (bottom panel). Dotted magenta line: maximum 
penetration of the convective envelope on the AGB (O.525M0). 
A dynamo operating on the RGB deposits magnetized material 
within mass shells below this line; above it, the magnetic field is 
sourced on the AGB. Reference time tio .9 ~ 12.37 Gyr subtracts 
MS and early giant branch evolution (R*(tio. 9 ) ~ IO. 9 R 0 ). 


of the remnant of an isolated, intermediate-mass star is 
therefore expected to grow and then gradually decay on 
this timescale (Section [6]) . 

3. MAGNETIC FIELD AMPLIFICATION NEAR 
CORE-ENVELOPE BOUNDARY 

We now investigate the growth of a magnetic field 
above and below the core-envelope boundary of a post- 
MS star. There are two profound differences here with a 
solar-type star. First, the radiative flux that is generated 
by the burning shell(s) much exceeds the flux emerging 
from the core on the MS. Second, hydrogen-rich material 
drifts downward rapidly through the convective bound¬ 
ary, as the core grows in mass. This leads to a secular 
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Fig. 2.— Interior profile of a Mzams = 5Mg star during RGB 
expansion and core He burning phase (top panel) and AGB phase 
(bottom panel). The deepest penetration of the convective enve¬ 
lope on the RGB remains outside the mass boundary of the WD 
remnant (~ 0.87 Mq). Dredge-up of helium occurs on the early 
AGB; the dotted magenta line shows the corresponding maximum 
penetration of the convection. 


accumulation of magnetic hclicity in the core, as this in¬ 
ward drift is combined with a persistent magnetic twist. 
Such a twist is generated when the magnetic field couples 
the rotation of core and envelope. In these respects, our 
considerations extend beyond existing dynamo theory. 


than about 0.1, corresponding to Co < 60. Specific ally in 
M-dwarfs, the closest MS analogs of giants. [Reiners et al.1 
(120091) observe that surface magnetic flux scales as Co 2 
for Co < 60. There does not appear to be a sharp cut¬ 
off of magnetic activity with Coriolis parameter, as the 
simplest mean-field dynamo models would suggest. 

The evidence therefore points to the presence of weak, 
but finite, dynamo activity at Co ~ 1 in MS stars. It also 
begs the question as to whether the relative strength of 
the dynamo-generated magnetic field (measured in terms 
of the convective stress at a given Coriolis parameter) 
differs in giants and MS stars. 

First it should be emphasized that a large-scale 
Maxwell torque that couples the rotation of the outer 
core to the envelope depends on a tiny poloidal field that 
is implanted from the dynamo layer. We quantify this 
minimal field in Section 13.31 

Some additional insight is provided by mean-field dy¬ 
namo theory, which suggests that the dynamo number 


D = 


allr 3 




(14) 


must exceed a threshold D > (27 r) 3 ~ 1 0 2 for magnetic 
field growth to be sustained (ICharbonneaulfMTS ). Here 
a is the parameter relating the toroidal magnetic field to 
the toroidal electromotive force, and v t is the turbulent 
diffusivity. The prospect for a dynamo is enhanced if the 
magnetic field can be anchored in a part of the star with 
relatively low (but finite) i/ t . 

This effect provides part of the motivation for consid¬ 
ering a thin, stably-stratified layer below the convective 
envelope as the locus of strong differentia l rotatio n, and 
a reservoir for a toroidal magnetic field (iParkerl I1993T) . 
However, the sharpness of the onset of stable stratifica¬ 
tion below the envelope boundary has presented a diffi¬ 
culty for this hypothesis. The solar convective motions 
have a relatively low Mach number (A4 CO n ~ 10 -4 ), and 
so the overshoot layer below the envelope is very thin. 
A toroidal magnetic field reaching equipartition with the 
convective stresses overcomes the stable stratification in 
a layer that is not much thicker. 

In Section 13.21 we show that the intense flux of radi¬ 
ation emanating from the core of a RGB or AGB star 
will induce an upward drift of magnetized material, one 
that more than counterbalances the downward flow to 
the burning shell(s). A tiny seed field that is generated 
by the convective motions will sustain this feedback (Sec- 
tion l3.3D . The interchange between toroidal and poloidal 
fields will persist in a layer of considerable thickness 
~ (0.01 — 0.1)^p over a dynamo cycle ~ (10 — 10 2 )r con . 
Convectively driven turbulence is strongly suppressed 
within this layer. 


3.1. Critical Co for Dynamo Action: 

Effect of an Intense Radiation Flux in the Tachocline 

Late-type MS stars generally are compact enough to 
sustain Co > 1 in their convective envelopes, as are sub¬ 
giants massive enough to have retained most of their na¬ 
tal angular momentum at the completion of core H burn¬ 
ing. Dynamo activity, as measured by chromospheric 
lines and coronal X-ray emission, begins to decline when 
the Rossby number Ro = P rot / r con = 27r/Co is larger 


3.2. Magnetic Buoyancy in the Tachocline 

We now examine the competition between the 
downward drift of hydrogen-rich material across the 
tachocline, and the upward drift of magnetized material 
(Figure [3]). Consider a horizontal (toroidal) magnetic 
flux tube, of a radial thickness comparable to the dis¬ 
tance Ar = -Rbenv — r to the base of the convective enve¬ 
lope. The flux tube starts in buoyancy equilibrium, with 
a density equal to that of its unmagnetized surroundings. 



























Fig. 3.— A magnetic flux tube positioned in the tachocline layer 
experiences an upward drift that is driven by the intense radiation 
flux through the tube. This drift can be strong enough to coun¬ 
terbalance the downward flow of hydrogen-rich material into the 
burning shell(s) below. 



r ( R o) 


Fig. 4.— Horizontal (toroidal) magnetic field anchored in a 
tachocline experiences radiatively driven buoyancy. Upward drift 
speed vf ?, Equation (ED, is compared with the downward flow 
of material from the convective envelope into the growing radia¬ 
tive core. Curves correspond to various timesteps near the tips 
of the RGB and AGB in our 1 Mq model, and one taken from 
the thermally pulsating AGB. The magnetic field is normalized to 
~ (47rpu^ on ) 1 / 2 , with vB proportional to B Buoyant drift ac¬ 
celerates toward the core-envelope boundary, due to the decrease in 
the Brunt-Vaisala frequency. Here downflow is evaluated directly 
from the stellar model. 


The temperature deficit in the flux tube (labelled ‘2?’) is 


Tb~T _ B\ 

T 8nP ’ 

which maintains a differential energy flux 

SF rad ^--^-(T B -T)T 3 

dup/Xr 


(15) 


(16) 


into the flux tube. This drives a buoyant rise of the flux 
tube at a speed v^, which is determined by 


R dS 

pk B TVr = 


dT rar ] d F r:l( \ 

dr A r 


(17) 


The background entropy gradient is usefully expressed in 
terms of the Brunt-Vaisala frequency 


so that 


N = 



7 — 1 dS 

7 dr ’ 



2^rad 

(A r) 2 pN 2 


(18) 


(19) 


The background flow of mass into the radiative core is 

dM core _ /| _ D 2 „/ d \ {dRcore „,nucA ^ 

dt ~ 4 ^benvP(^be„v) ( ^ J ~ £ . • 

( 20 ) 

After the star passes onto the thermally pulsating AGB, 
the relation between core growth rate and luminosity L 
only holds after averaging over pulsations. The nuclear 
energy released per unit mass of new core material is 
E\ = 6.3A MeV/m„ during the first step of converting 
material with hydrogen mass fraction X to helium; and 
£2 = Si + 0.8 MeV/m„ during double shell burning. The 
corresponding velocity is 


v 


nuc 

r 


Arad 

p£i 


( 21 ) 


The ratio vf /\v^ uc \ is plotted in Figure 0] at various in¬ 
tervals during the RGB and AGB evolution of our 1 Mq 
model. One observes that the most strongly magnetized 
parts of the tachocline will not be subducted if the mag¬ 
netic pressure approaches the kinetic pressure at the base 
of the convective envelope. It should be kept in mind 
that this effect depends on a difference in magnetization: 
subduction of the less weakly magnetized parts of the 
tachocline must continue even where v^ > [n" uc |. 

Magnetic buoyancy has the effect of expanding the ac¬ 
tive dynamo region below i?benv, by allowing a layer with 
strong dCl/dr - the tachocline - to remain in contact with 
the convection zone0 We can estimate the size of this re¬ 
gion by comparing the buoyant rise time A r/v^ (r) with 
the dynamo timescale -Pdyn (set equal to 100r con ) 


Ar = Abenv - r = v^(r) ■ P dyn . (22) 

To invert this expression, we fit (Ar) to a power law 
in Ar, measuring it in two different radial zones. 

Figure [5] shows the resulting buoyant shell thickness, in 
comparison with the local scale height, during the AGB 
phase of the 5 Mq model. The buoyant shell thickness 
reaches ~ 0.1Zp(i?benv) before and during the first ther¬ 
mal pulse, and then relaxes to a value about 10 times 
smaller. 


3.3. Seeding a Toroidal Magnetic Field 
in the Tachocline 

We now consider the minimal seed field that will trigger 
dynamo feedback in the tachocline, thereby facilitating 
the magnetic coupling of core and envelope. We suppose 

2 We do not address here whether the re-conversion of wound 
toroidal magnetic field back to the poloidal direction is concen¬ 
trated in the tachocline, or instead the convection zone, as in 
Parker’s 1993 model. 
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t (years) 

Fig. 5.— Thickness A r of magnetized shell below the core¬ 
envelope boundary which remains in contact with the inner con¬ 
vection zone over the dynamo period Pdyn- See Equation H22D . 
Here A r is normalized to the local pressure scaleheight. Plot cov¬ 
ers the thermally pulsating AGB phase of the 5 Mq model, with 

B <t> = ( 47r W?on) 1/2 - 

that this seed field is generated within the convective 
envelope, and is deposited in the outer radiative core by 
the mass flow toward the burning shell(s). 

This seed field must pass a certain threshold if lin¬ 
ear winding in the tachocline is to generate a buoy¬ 
antly unstable toroidal field. The time available for 
linear winding depends on the depth Ar below the 
radiative-convective boundary. Consider an unmagne¬ 
tized core in which matter flowing across the convec¬ 
tive boundary sustains an angular velocity profile 17 = 
17benv(?V-ffbenv) -2 - Near the boundary, the toroidal mag¬ 
netic field that is sourced by a seed radial field B see <j 
is B ^ = —2B see dl7benvAr/u™ c . Substituting this into 
Equation m and setting > |u/ uc | gives 

BL d N- P\vr c | 3 

47r/3(i?benv)'C?on 2 M 2 V% on F rad ' 

This minimal seed field for dynamo feedback 
turns out to be quite weak. Figure [6] shows 
B S eed/s/ 47rp(i?benv)'c C on during the final AGB evolu¬ 
tion of a 5 Mq star, measured at a distance Ar = 
/p(i?benv)/10 below the core-envelope boundary. Much 
of the C/O mass deposition from double shell burning 
is concentrated around t ~ 1.0716 x 10 8 yr, when this 
quantity is as small as ~ 10" 9 ; it rises to ~ 10 -4 during 
the phase of rapid thermal pulses. 

Hydrodynamic turbulence is strongly suppressed in the 
tachocline. The mean-field approach to dynamo action 
then suggests (see the review in Section m that dy¬ 
namo activity should be easily sustained in the tachocline 
as long as the toroidal magnetic field can be rotated 
back into the poloidal direction. This effec t is present 
in numerical simulations of a shear layer dCline et al.1 
I2003T) : and is conjectured to operate in a more global dy¬ 
namo model where toroidal magnetic field formed in the 
tachocline diffuses into the convective envelope (1 Parked 



t (years) 

Fig. 6. — Minimal seed poloidal magnetic field, deposited in the 
outer core of an AGB star, which will grow by linear winding 
to a sufficient strength that radiatively driven magnetic buoyancy 
overcomes the dow nwa rd drift of envelope material to the burning 
shell(s). Equation ll23l is evaluated at depth Ar = Zp ) benv/10 be¬ 
low core-envelope boundary. Plot covers same thermally pulsating 
AGB phase of 5 M® model as Figure (5) 



Fig. 7. — The magnetic helicity (expressed as a magnetic flux 
AH 1 / 2 ) that accumulates in the time tdrift — £p/ |u£ uc |. RGB 
(left) and AGB (right) evolution of Mzams — 1 Mq model with 
a Jupiter companion in an initial orbit a,i = 1 AU. Most of 
the helicity generated on the RGB is concentrated near the tip 
(orange short-dash curve shows cumulative helicity). This gives 
Mbenv(5.65 x 10 7 ) ~ 0.388 Mq as the effective base of the magne¬ 
tized region. The net mass passing through the dynamo-active 
layer is then ~ O.12M0 on the RGB, and ~ 0.014M® on the 
AGB. The grea ter helicity accumulated through the ‘f2(0)’ channel 
(Eq uation l‘27l) h as compared with the ‘inflow’ channel (Equation 
J33J), is reflected in the final WD magnetic field (Figure [9]). 


GUI. Estimates of the equipartition toroidal magnetic 
field and the large-scale Maxwell stress are given by equa¬ 
tions CH and m ■ 
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Fig. 8. — Comparison of magnetic helicity generated over thrift = 
£p/\v^ uc \ through the same two channels as Figure 0 but now in 
the Mzams = 5 Mq model without a planetary companion. Only 
the AGB contributes directly to final WD magnetization in this 
case, within a mass shell ~ 0.044 Mq. The largest net contribution 
comes from the first thermal pulse. 


4. WHITE DWARF MAGNETIC FIELD 

After magnetic helicity accumulates in the core, we 
allow the magnetic field to relax to a more isotropic con¬ 
figuration. A first estimate of the flux threading one 
hemisphere is 

$ core ~ H 1/2 . (24) 

The final poloidal magnetic field of a WD of radius R w d 
is 

5wd ~ (25) 

We first evaluate the magnetic helicity that is de¬ 
posited in a post-AGB star, using the 1 Mq and 5 Mq 
stellar models constructed using MESA. This stored he¬ 
licity is then converted to a dipolar magnetic field in the 
WD remnant using Equations (|24l) and (|25|). 

This estimate presupposes that toroidal and poloidal 
magnetic fluxes evolve to an energy minimum with <f> p ~ 
$0, all the while preserving the product helicity TL ~ 
Such an interchange may occur later in the life 
of a WD as the result of ohmic diffusion (Section 0; or 
earlier on as the result of relaxation within episodes of 
confined convection (such as the core helium flash, or 
helium shell flashes). 

4.1. Accumulation of Magnetic Helicity in a 
Radiative Core Relaxing to Solid Rotation 

The zone of interest here is the tenuous mantle of radia¬ 
tive material that is sandwiched between the hydrogen- 
depleted core and the base of the convective envelope. 
Here only a small mass < 10~ 3 Mq lies within a scale- 
height, and so we model the mass inflow through the 
mantle as quasi-steady, 

M cole = 47rr 2 |u r |p(r). (26) 


The poloidal magnetic field is approximated as ra¬ 
dial with flux function 0, and the toroidal magnetic 
field is expressed in terms of a radial magnetic twist as 
Bj,(r,8) = {d<t>B/dr)r s\\\9B r . Then helicity accumu¬ 
lates in the core at the rate 


dt 


—2tt Sv r 


d9$ r 


d$ r d(j) B 
89 dr ’ 


(27) 


where 5v r is the flow velocity relative to the core-envelope 
boundary. 

Now consider the torque imparted to a mass shell by 
the Maxwell stress B r B^/4 -7T, 


— = S 

dt\ 3 J 


dS r sin 9 


B r B 4 j 
47T 


(28) 


Here the Lagrangian time derivative follows the radial 
flow at speed v r - Since d(Sm)/dt = 0, one has, following 
Equation 0, 


d 

dt 


- cn = 


87 ipr 2 


d9 sin 9 


(9<E> r 

~89 


d 2 (j)B 

dr 2 


(29) 


Both hemispheres generally contribute in the same sense 
to the torque. 


4.1.1. Helicity Sourced by Radial Inflow 

We suppose that the magnetic field is strong enough 
that nearly solid rotation is maintained in the radiative 
layers, \dCl/dr\ <C f l/r. Each mass shell starts with a 
mean angular frequency flbenv that is imposed by the 
lower part of the convective envelope. The twist that is 
required to maintain a weak angular velocity gradient is 
found from Equation (l29l) . 

8 f>B 2Af core ribenv^ 

~87 2 (. B r r 2 ) 2 ’ 1 j 


where the coefficient on the right-hand side corresponds 
to B r independent of 9. From Equation (1A10I) one has 


d(j)B r 8Q 
dr | | dr 


(31) 


Nearly solid rotation is maintained if 

d 2 

2 D r ^ 2 

s 3^ >> 


(32) 


Now consider the magnetic helicity that is advected 
into the growing radiative core. We evaluate the mag¬ 
netic twist in Equation (1271) from Equation (l30l) , obtain¬ 
ing 


— = /(AM core ) • 7rM core n b env-Rbenv^W(-Rbenv) (33) 

at r — 7?benv 

The factor /, with indeterminate sign but magnitude 
|/| < 1, represents the cancellation between hemispheres. 
We discuss its origin in Section 14.1.31 


4.1.2. Magnetic Twist Compensating a Latitude-dependent 
Convective Torque 

A greater Maxwell stress results from the application 
of convective stresses to the outer radiative core. A con¬ 
vection zone generally sustains latitudinal gradients in 
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rotation. Part of the differential Reynolds stress is trans¬ 
mitted through the tacliocline to the core by a dynamo¬ 
generated magnetic field (Equation (1121) 1. As a working 
estimate, we use 


—- = £sp(Rbenv)(^benv-Rbenv) 2 , (34) 

47T 

with £b ~ 10 3 , corresponding to a dynamo period 
Pdyn ~ 10[Afl(l?benv)] _1 and a pole-equator offset 
Af2(77benv) ^ f^benv- 

Given that the deeper parts of the radiative core are 
coupled magnetically to the tacliocline, uniform rotation 
can be maintained in the outer core only if this stress is 
divergence-free below the tacliocline. The corresponding 
magnetic helicity flow into the core through the advection 
of this twisted field is 


<m 

dt 


/(AM core ) • TT£ B M coie R h env (fi b env Rh env ) 2 . (35) 


The contributions to 77 from equations (1351) and (l35l) 
are labelled in the Figures as ‘inflow’ and ‘f1(0)’, respec¬ 
tively. 


4.1.3. Cancellation Between Hemisphere 

One feature of Equations (l33l) and (l35l) stands out: the 
growth of H is independent of the seed radial magnetic 
field, as long as it is strong enough to satisfy the inequal¬ 
ity (1321) . This means that the contribution to A ■ B in one 
hemisphere is not strongly affected by fluctuations in the 
magnitude and sign of the poloidal magnetic field that 
is supplied to the radiative layer by a dynamo operating 
near the convective boundary. 

A net mass AM core flows through the convective- 
radiative boundary over the duration of dynamo activity. 
In the case of an AGB star, this is approximately the in¬ 
crease in the mass of the C/O core due to double-shell 
burning. Note that AM core is generally much larger than 
the mass of the rarefied core mantle (~ 1O~ 3 M 0 ). The 
low mass of the mantle imposes a short-timescale cutoff 
to fluctuations in the summed 77 from the two hemi¬ 
spheres. We express this cutoff in terms of the mass 
within a scale height below the core-envelope boundary, 


AM mi „ 


47rPr 4 


GM C , 


Rbe 


Then we take 


/ AM • \ 1//2 

mT 


(36) 


(37) 


4.2. Results 

First consider the M z ams = 1 M 0 model. The an¬ 
gular momentum profile of the star is evolved in post¬ 
processing using the method of Section [Til and Paper I. 
Companion planets with a range of masses (Earth, Nep¬ 
tune, and Jupiter) are placed at an initial semi-major 
axis a* = [0.5,0.75,1,1.5, 2] AU and with a range of or¬ 
bital eccentricities. 

Magnetic fields generated during the RGB and AGB 
phases both contribute to the remnant WD field, but at 
different depths in the star (see Figure [T]). A net mass 
~ O.135M 0 passes through the core-envelope boundary 


while conditions are favorable for magnetic field growth. 

0 

Figure [T] shows the helicity that accumulates in the 
brief period idrift = ^p/|u" uc | when stellar material set¬ 
tles through the outer scale height of the radiative core. 
The dominant contribution comes from the ‘f2(0)’ chan¬ 
nel. 

In the 5M 0 star, only helicity accumulated on the AGB 
contributes directly to the final WD magnetic field (Fig¬ 
ure 0. During this phase ~ O.O44M 0 passes through the 
active dynamo region (given our normalization r]B = 0.05 
of the Blocker mass loss formula). We see in Figure[5]that 
the initial phase of the dynamo, which occurs during the 
first thermal pulse, contributes the majority of 77. In the 
5 Mq model, we simplified the treatment of the thermal 
pulses by linearly interpolating all of the relevant quan¬ 
tities between the local minima in 7?benv after each He 
shell flash, so as to focus on the mean growth of core 
mass and 77. The stellar model considered in Figure [5] 
has no planetary companion. 

The resulting WD dipolar field is shown in Figure [9] 
for both the 1 an d 5 M m m o dels, usin g the mass- 
radius relation of lZapolskv fc Salpeter (196 9). A strong 
magnetic field in the 0.55 M 0 WD remnant of the 1 M 0 
progenitor depends on the injection of angular momen¬ 
tum from a planetary companion (as considered here), 
or a tidal interaction with a stellar companion. We find 
that 77 wc j is an increasing function of the angular momen¬ 
tum absorbed; the various points correspond to different 
realizations of the orbital evolution. 

Dipole fields approaching 10 s G are achievable through 
the ‘£2(0)’ channel, following the absorption of a Jupiter- 
mass planet. Note that an Earth-mass planet does not 
signicantly augment the angular momentum left over at 
the end of the MS spindown, and so this case corresponds 
closely to an isolated star that does not interact with 
planets. 

Figure 0 focuses on those runs which end in the engulf- 
ment of the planet. However there are cases in which the 
planet is pushed out into a large orbit due to the con¬ 
vective quadrupole. The planet therefore gains angular 
momentum from the star, which in turn causes the star 
to spin in the opposite direction. In some cases this leads 
to the development of a dynamo as well. We ignore these 
particular cases. 

An external source of angular momentum injection is 
not required to sustain a dynamo in the 5 M 0 progenitor. 
Nonetheless, B w d depends on the C/O mass that accu¬ 
mulates during double shell burning, and therefore on 
the mass-loss rate during the expulsion of the envelope. 

The effect of adding a relatively massive companion to 
the 5 Mq model is considered in Figure HUl The absorp¬ 
tion of a 3 Mj planet significantly augments the remnant 
dipole field in the 0.87 M 0 WD. We see a steady increase 
in the remnant field as the companion mass increases to 
the brown dwarf range, or even to ~ 0.1 M 0 , for the 
angular velocity profile considered here (a = 3/2). For 
comparison we show the magnetic fields obtained by a 
shallower inner rotation profile (a = 1), which may be 
applicable depending on the strength of the Coriolis back 
reaction. In this case an increase in angular momentum 

3 For example, we do not include the first dredge-up phase, dur¬ 
ing which the convective envelope grows in mass. 
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Fig. 9.— Dipole magnetic field, Equation m- left behind in the WD remnant of a 1 Mg star interacting with a planet (black = Earth 
mass, red = Neptune mass, blue = Jupiter mass). Result is plotted separately for the two helicity channels described in the text. Horizontal 
axis shows initial orbital angular momentum of planet companion. Each tick represents an average over realizations of a given set of orbital 
initial conditions (a; and e,). Magenta squares show result for 5 Mg progenitor with initial equatorial rotation speed 50 km s~ 1 and no 
companion; orange squares same model with angular momentum of 3 Mj, a, = 2 AU planet added during early AGB. Dynamo shuts off 
when Co benv < Co benViCrit = 0.1 (left panel) or Co benVjCrit = 0.3 (right panel). 



Fig. 10. — Dipole magnetic field obtained via the £1(9) channel 
in the 0.87 Mg remnant of a 5 Mg star that absorbs an orbiting 
companion (a; = 2AU) of various masses during the early AGB 
phase f R r = 200i?g). Top points: rotation profile £i(r) oc r —1,5 

S -2 ) in the inner (outer) convective envelope (« = 1.5 in Equation 
). Bottom points: a = 1. 


results in a decrease of remnant magnetic field, caused by 
a decrease in Cobenv due to the shallow T CO nv(?') profile. 

5. WHITE DWARF ROTATION 

We now analyze the spin evolution of the core and 
envelope of our stellar models. Our focus is on the mag¬ 
netic de-coupling between the core and envelope. Three 
factors determine the effectiveness with which a poloidal 


magnetic field transfers angular momentum across the 
core-envelope boundary, and establishes solid rotation in 
the outer core: i) the timescale over which mass is ex¬ 
changed between core and envelope; ii) the direction of 
this exchange; and iii) the Coriolis parameter in the inner 
envelope. 

The internal structure of the star evolves slowly enough 
during the early ascent of the giant branches, and during 
core helium burning, that a magnetic coupling between 
core and envelope is very difficult to avoid (see Figure 2 
of Paper I). There is a much more rapid exchange of mass 
between envelope and core near the tips of the RGB and 
AGB. The first giant phase is (except in cases of close 
binary interaction) followed by more extended evolution. 
So our focus here is on the final AGB superwind phase 
during which most of the envelope mass is removed. 

The ability of the inner envelope to generate a mag¬ 
netic field is most relevant while the radiative core is 
growing in mass. Then magnetized material is advected 
downward into the outer core, where it can communi¬ 
cate changes in rotation rate in the envelope to the inner 
core. The core grows in mass during the thermally pul¬ 
sating AGB phase, after averaging over the pulsations; 
but following the contraction of the envelope the super¬ 
wind causes a small but rotationally significant decrease 
in mass. 

An important effect involves the re-expansion of the 
envelope following a late thermal pulse: following this 
about 10 -3 Mq of hydrogen-rich material is transfered 
to the convective envelope. If core and envelope are still 
magnetically coupled, this causes a dramatic spindown 
of the core through the internal exchange of angular mo¬ 
mentum. 

Faster WD spin is mediated by an early loss of angular 
momentum from the surface of an AGB star. This loss, 
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Fig. 11.— Rotational angular momentum of Mzams — 1M© 
AGB star (dotted red line) and its core (short-dash orange line) just 
before and during the super-wind phase. The star has absorbed 
a Jupiter-mass companion near the tip of the RGB. Black line: 
Coriolis parameter at base of convective envelope. Rotation profile 
fl(r) oc r ~ 1,5 (r~ 2 ) in the inner (outer) convective envelope. When 
calculating remnant WD spins, we freeze Q cor e when and if Cob env 
drops below a critical value, beyond which time the curves plotted 
here do not apply. 



if strong enough, will push the Coriolis parameter Cobenv 
of the inner envelope below unity well before the envelope 
contracts. The remnant WD spin period depends on the 
critical value of Cobenv below which the dynamo effec¬ 
tively Shuts Off. If the transition tO CObenv < CObenv, crit 
takes place before the final shell flash, then the angular 
momentum of the core is effectively frozen in and inher¬ 
ited by the WD. On the other hand, if this transition 
occurs afterward, then we expect that the rotation of 
the outer core remains coupled to the inner envelope as 
the photosphere of the star contracts. 

Figures fill and fl2l show the evolution of the net stellar 
angular momentum J* and the core angular momentum 
J C ore during the AGB phase, under the assumption of a 
tight core-envelope coupling. The internal rotation pro¬ 
file, the prescription for mass and angular momentum 
loss, and the interaction with a planetary companion, 
are the same as those described in Sections iii&izi and 
detailed in Paper I. As the star experiences strong mass 
loss, there is a rapid drop in J*. 

During this phase, J* is still dominated by the en¬ 
velope. Also plotted is Cobenv for comparison. Core¬ 
envelope decoupling becomes plausible at Cobenv = 
Cobenv,crit ~ 0.1-0.3, given the strong dependence of 
magnetic field on rotati on that is obse rved in deeply con¬ 
vective MS stars dReiners et al.ll2009l) . 

Whether the transition to Cobenv < Cobenv,crit is en¬ 
countered before the final He shell flash depends on both 
J* and the value of Cobenv,crit- The transition occurs 
earlier for a larger critical Coriolis parameter, as in our 
1 Mq model with Cobenv = 0.3. But we find a tighter 
restriction on J* when Cobenv is reduced to 0.1: then the 
absorption of a Jupiter will cause core and envelope to 
remain coupled past the final flash in the 1 Mq model, 
but the absorption of a Neptune (or of no planet) al¬ 
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Fig. 12.— Similar to Figure fTTl but for the Mzams = 5 Mq star. 


lows a transition to a weakly magnetized envelope. We 
find that Cobenv generally remains above Cobenv,crit at 
the final shell flash in the 5 Mq model, implying a tight 
core-envelope coupling during the transition to a post- 
AGB star. 

In the following two sections, we separately consider 
the rotational evolution of our stellar models with, and 
without, a transition to weak core-envelope coupling. 

5.1. Dependence of White Dwarf Spin on 
Coriolis Parameter for Core-Envelope Decoupling 

The WD spin period that results from core-envelope 
decoupling at Cobenv crit = 0.1 and 0.3 in the 1 Mq model 
is shown in Figure fTHl For the more conservative estimate 
Cobenv,crit = 0.1, we find that P wc j ranges over 0.5-1.5 d. 

In this situation, where the timing of the transition 
to weak magnetization in the envelope depends on J*, 
P w d depends on the initial orbital angular momentum of 
the planetary companion. A lower total angular momen¬ 
tum results in faster decoupling, and therefore a larger 
trapped core angular momentum: the core rotation fre¬ 
quency at decoupling is tied to the convective timescale 
in the inner envelope, and to the structure of the star. 

This behavior is explained in Figure 1141 which shows 
the minimum angular momentum required to maintain 
Cobenv = 0.1 in the 1 Mq model near the tip of the AGB. 
This is compared with the actual angular momentum 
evolution that results from the absorption of a Jupiter, 
Neptune or Earth, initially orbiting at a* = 1 AU with 
eccentricity e* = 0. 

5.2. White Dwarf Spin in Case of Continued 
Core-Envelope Coupling 

We observe a significant difference in the remnant an¬ 
gular momentum in two models which maintain a tight 
core-envelope coupling through to the post-AGB phase: 
the 5 Mq model, and the 1 Mq model which absorbs a 
Jupiter-mass planet. In the first case, downward angular 
momentum pumping in the convective envelope allows 
</ CO re to remain relatively constant during the ejection of 
the envelope, and only modest spindown is caused by a 
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Fig. 13.— Rotation period of 0.55 Mq WD remnant of 1 Mq star, 
as a function of the initial orbital angular momentum of planetary 
companion. Squares: Core and envelope remain magnetically cou¬ 
pled until the post-AGB phase, x’s: Core-envelope decoupling be¬ 
fore the final helium shell flash, with J w u set to J co re at decoupling. 
Decoupling takes place when Cobenv drops below Cobenv crit = 0.1 
(upper panel) or 0.3 (lower panel). Horizontal magenta line: re¬ 
sult for 0.87 Mq remnant of Mzams — 5 Mq star with 50 km s —1 
equatorial rotation period. 
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wind from the surface of the post-AGB star (Figure [12]). 
The spin of the WD remnant is directly proportional to 
the initial spin period of the star, except in the case of the 
absorption of a relatively massive (> 3 Mj) companion. 
An equatorial rotation speed 50 km s -1 on the zero-age 
MS results in P w( j ~ 0.6 days. 

The initial contraction of the envelope also does not 
cause much of a reduction in J core in the 1 Mq model; 
but we observe a much stronger decrease following the 
re-expansion and re-contraction of the envelope that is 
triggered by a late helium shell instability ('Figure fill) . 
This means that the WD remnant of the 5 Mq star has 
a spin period P w d ~ 15 hr, as compared with ~ 2 yr for 



Fig. 14.— Black solid line: minimum J* that maintains Cobenv = 
0.1 near the tip of the AGB of a Mzams = 1 Mq star. For com¬ 
parison, colored lines show calculated angular momentum of star 
absorbing a Jupiter (red dotted), Neptune (orange short-dashed) 
or Earth (blue long-dashed) starting with a; = 1 AU, e,; = 0. J CO re 
is shown separately for Neptune and Earth interaction, where core 
and envelope decouple relatively early. Decoupling is delayed to 
the final super-wind phase when the companion has a Jupiter m ass, 
resulting in lower J co re and the longer WD spin seen in Figure [13] 


the remnant of the 1 Mq star. 

To understand how such a dramatic difference in fi¬ 
nal WD spin could arise, we consider a simplified ana¬ 
lytic model of mass and angular momentum loss from an 
AGB star. The envelope of the model star has uniform 
specific angular momentum, the core has a fixed radius 
-Rbenv, and the core rotates as a solid body with angu¬ 
lar velocity f2 core = Ubenv- Mass is lost from the outer 
boundary of the star, and is also exchanged between core 
and envelope, 


dM env _ dM+ _ dM henv 
dt dt dt 

Here M env is the mass of the convective H-rich envelope, 
and Mbenv is the mass of (mainly radiative) material in¬ 
side the base of the envelope. The moment of inertia of 
the core is initially a small fraction of the total effective 
moment of inertia 


Pff 


T 2 

P 2 _i_ T 

q A'knTnbenv t j ci 

^ ^benv 


(39) 


where we parameterize I COTe = (2e cole /3)M henv Rl env . 

The coefficient e cor e is very small, but a key consid¬ 
eration is its relative size compared with the fraction 
M env /M* of the stellar mass that is contained in the en¬ 
velope. We find that e cor e ranges over ~ 3 x 10” 5 — 10 -3 
during the superwind and post-AGB phases. The en¬ 
velope of the 1 Mq model completes its first contrac¬ 
tion while it retains ~ 10~ 3 Mq, but following its re¬ 
expansion and re-contraction this mass has dropped to 
M env ~ 1O _4 M0. In the 5 Mq model, the envelope ex¬ 
periences a single contraction when its mass drops below 
~ 10 _4 M q . See Figure [T5l 

Consider first the stage(s) when the envelope is still 
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inflated, and the core radius i?benv changes only slowly. 
Then 


and 


dJ* 2 q . cLM* 


__ d2 q 

0 -^benv^^benv" 


dMf, 


dt 


A-f ■ 




benv 


dt 


1 dfl 


benv 


dM benv/ dt 


dt 




(40) 

(41) 


Changes in the radius of the star do not enter into these 
expressions as a result of the uniform specific angular 
momentum profile in the envelope. 

The core mass increases relatively slowly compared 
with the rapid drop in M env during the peak of the super¬ 
wind. So flbenv increases only slightly up to the contrac¬ 
tion of the envelope, where it takes the value flbenv.coi- 
Beyond that point, the entire star rotates nearly as a 
solid body. A wind from its surface, which carries away 
a mass AM* causes a net spindown 


^benv,col 
1-benv.col 


~ exp 


1 AM* 

Ccore A/* 


(42) 


This works out to a factor ~ 0.2 decrease in fibenv in the 
1 Mq model at the onset of the late helium shell flash. 

The key step in the dramatic spindown of the post- 
AGB star involves the re-expansion of the envelope. This 
occurs by the transfer of ~ 10 -3 M 0 of hydrogen-rich 
material from the radiative part of the star into a reju¬ 
venated convective layer, all occuring at nearly constant 
total angular momentum. Now the core rotation rate 
(still equal to the rotation rate at the base of the enve¬ 
lope due to a tight magnetic coupling) decreases to 
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fi? i M el 

benv,col ei 


e Mbenv 


eMbe 
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30' 

(43) 


(Although the star contracts by a factor ~ 10 in between 
the initial envelope contraction and the late helium shell 
flash, e CO re returns to its pre-contraction value after the 
re-expansion of the envelope.) The rotation rate of the 
core remains nearly constant at the value (1431) . as most 
of the rejuvenated envelope is expelled. 

The net spindown from the initial contraction of the 
envelope works out to flbenv/^benv.coi ~ 1/150. This 
behavior is neatly encapsulated by the trend of J* and 
Jcore with the mass M* — M w d remaining to be expelled 
(Figure dS]). 


5.3. Minimal Magnetic Flux Enforcing 
Strong Core-Envelope Coupling 

Here we evaluate the minimal magnetic flux that must 
thread the core-envelope boundary in order for i) the 
outer radiative core to remain in solid rotation; and 
ii) the rotation of the inner core to remain coupled 
to the outer core and inner envelope. A related (al¬ 
though less quantitative ) estimate has been made by 
iSnruit fe Phinnevi (119981 ) for the central cores of massive 
stars during the later stages of nuclear burning. 

First consider the outer core. Although it contains a 
small fraction of the core mass, it can contribute signif- 



M star -M wd (M 0 ) 



Fig. 15. — Total spin angular momentum of AGB star during 
the contraction of its hydrogen envelope (black solid line) and the 
angular momentum of the core (red dotted line). Top panel: 1 Mq 
model, including both the initial envelope contraction and the re¬ 
expansion following a late helium shell flash. Bottome panel: 5 Mq 
model. We show for comparison the dimensionless core moment of 
inertia e CO re = icore/(2/3 )M benv l?g env (orange short-dash line). 


icantly to the moment of inertia. We define a charac¬ 
teristic magnetic flux threading a single hemisphere by 
Setting B r / \J 47rp(i?benv) = I Vr dRb env /dt\'. 

^solid = ^"-^benvl^r di^benv/ dt \ \J 47Tp(i?benv) • (44) 

Next consider changes in the angular momentum of the 
entire core, which are sensitive to the evolving rotation 
profile of the envelope. The core angular momentum re¬ 
sponds to the applied Maxwell stress B r B^I 4-7T according 
to 

Icoi e d ^° IC ~ J Rbenv sin 6 ■ 2tt Rl env sin Odd. 

( 45 ) 

Taking a second time derivative, making use of the in- 
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duction equation, we get 


d 2 fl r 


dt 2 


- / 


dr 


, sin 3 Odd. 


(46) 


Here the poloidal magnetic field is assumed to be purely 
radial, so that d r (r 2 B r ) = 0, and the integral is evaluated 
at radius Rbenv Integrating over 6 then gives 


d 2 fL 


dt 2 


lBlR 4 


an 

benv ~dr 


(47) 


We take the rotation rate to vary according to 


^(r, t) ~ f^benv 



(48) 


Then to avoid spindown on a timescale tq , the magnetic 
flux must exceed 
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We note that tq is controlled by the rate at which mate¬ 
rial is ejected from the envelope during the final super¬ 
wind phase, and as a result i’coupie can in principle be 
larger or smaller than ^soiid (which is determined by the 
requirement that the core itself remain in solid rotation). 

These estimates of the magnetic flux that will enforce 
a tight rotational coupling between core and envelope 
are compared in Figures [TB] and [T7] with the radial flux 
(AH) 1 / 2 that accumulates at the core-envelope bound¬ 
ary. This increment AH is calculated using Equation 
(l35l) . taking the timescale for helicity accumulation to be 
the minimum of the radial drift time £p/ |i;“ uc | and tq. 

The larger value of the dynamo-generated flux sup¬ 
ports our choice of Cobenv < 1 for c ore-envelope decou¬ 
pling. Recall also from Section 13.31 that a tiny poloidal 
magnetic field - corresponding to 10~ 4 y / 47r/9(Hbenv)'Ccon 
during the later stages of the AGB will seed a toroidal 
field in the tachoclinc that is strong enough to interact 
buoyantly with the layers above (see Figure [6]). 


6. MAGNETIC FIELD EMERGENCE AND DECAY 

We now consider the emergence of magnetic fields from 
post-AGB stars, and their subsequent decay. The life¬ 
time of the visible magnetic field at the surface of the 
star depends on the thickness of the magnetized layer, 
which in turn depends on the mass of the progenitor and 
the mass loss rate during the expulsion of the hydrogen 
envelope. 

The strong dredge-up of helium that is experienced by 
stars of mass >2.3 Mq leaves only a thin mass shell to be 
processed on the AGB. The remainder of the C/O mate¬ 
rial is generated during core He burning, when the central 
convective core is decoupled from the exterior and there¬ 
fore must conserve H. On the other hand, the helium 
generated on the RGB is not dredged up in solar-mass 
stars, and so magnetic fields deposited during both the 
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Fig. 16.— Comparison between magnetic flux r T’ C o Up ]e (Equation 
1491 . red dotted line) that mai ntai ns core-envelope coupling on 
the AGB, flux < I > so iid (Equation 1 44 ll . purple dot-dashed line) that 
maintains solid rotation in the outer core against the inward mass 
flow, and flux AH 1//2 generated by the 0(0) dynamo process on 
radial drift time (orange short-dashed line). Blue long-dashed line: 
magnetic flux corresponding to a poloidal field in equipartition with 

the convective motions, <4>eq = rri? 2 env V 4?r P(^benv ) v con- Solid 
black line: Coriolis parameter at base of envelope. We argue that 
the strong dependence of the dynamo-generated field on Cobenv 
allows for core-envelope decoupling below Cobenv ~ 0.1 — 0.3. 
Fluxes are only plotted up to the final thermal pulse: afterward 
dM\> env /dt < 0 and the coupling of the outer core to the envelope 
is less sensitive to the instantaneously generated magnetic field. 
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Fig. 17. — Similar to Figure ITBl but for the Mzams = 5Mq 
model. 


red giant and asymptotic giant phases contribute directly 
to the remnant WD field. 

We found that mass loss on the AGB can push the con¬ 
vective envelope below the threshold for dynamo activity 
(Figures fTTI and [T^]) . In this case, a thin outer layer of 
the remnant is initially unmagnetized; a strong magnetic 
field only emerges by ohmic diffusion. We do not consider 
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any contribution to the magnetic field from the contrac¬ 
tion of the final 10 -4 M 0 of envelope material, which 
might be rapidly rota ting as the result of anisotropic 
mass loss (jSprnitl 11998h . 

The timescale for flux emergence, and the asymptotic 
decrease in the trapped magnetic flux, are easily obtained 
from the induction equation 0 


dB 

~dt 


— V x ( 77 V x B). 


(50) 


Here 77 = c 2 / 47 rcr is the magnetic diffusivity expressed in 
terms of electrical conductivity er and speed of light c. For 
degenerate matter in the liquid state the conductivity can 
be es timated using (lYakovlev fc Uroinl 119801 : lltoh et ahl 
1983), 


a = 8.5 x 10 21 s _1 


1 x 3 
A ei (Z) 1 + x 2 ' 


(51) 


Here A e j ~ 1 is the Coulomb logarithm; (Z) = 
XiZf /Ai, with Xi the mass fraction, Zi and Ai 
the nuclear charge and mass of species i, and /i e is the 
mean molecular weight per electron. Finally x = pf/T n e c 
measures the Fermi momentum of the electrons (of mass 
m e ). A good approximation to the diffusion timescale 
through a pressure scale height Ip is 


^diff ,/ p 


Analp 


(52) 


An analytic approximation to fdiff.ip is easily obtained 
in a thin outer shell of mass AM that is supported by 
the pressure of non-relativistic electrons, 


P = Kp 5/3 
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(54) 


in a nearly uniform surface gravitational field g. The in¬ 
tegration of the equation of hydrostatic equilibrium gives 
a pressure scale height Ip = (2/5 )(i? w d — r) and density 
p(lp) = ( glp/K ) 3//2 . Substituting a from equation m 
leads to 


10 in io f AM\ 7/5 f M wd \~ 3/5 
idiff,analytic ~ 1.8 x 10 years y M@ J [mq) 

flwd y S/5 (YA 2 {(AY 1 A -1 
io- 2 i? 0 ; \0.5y { 6 ) ei ■ 

(55) 


From this expression it is easy to see that the diffusion 
time through a shell of fixed AM depends weakly on the 
total WD mass. 

The diffusion time through the magnetized surface 
layer of our 0.55 and 0.87 M 0 model WDs is shown in 
Figures [18] and [T9] The longer diffusion time in the 
0.55 Mq remnant is due to the larger magnetized mass. 

A direct integration of the induction equation is 
straightforward if we ignore the continuing hydromag- 
netic adjustment that must accompany ohmic diffusion. 


4 Hall drift does not modify an axisymmetric magnetic field in 
a fluid star, and so we consider only ohmic effects here. 



Fig. 18.— Ohmic diffusion time in magnetized l ayer s of 0.55 Mq 
WD remnant of 1 A /.., s tar, as given by Equations (1521) (black line) 
and the analytic fit 1551 (magenta dashed line). Mass shells passing 
through the dynamo-active layer on the RGB (AGB) are bounded 
by green (blue) dashed lines (when the star has absorbed a Jupiter- 
mass planet). The hump in the diffusion times close to the surface 
is caused by a shift in composition and electrical conductivity. 


Making the same approximation of a geometrically thin 
magnetized layer, one has 

d®r,<t> _ c 2 d 2 $ r , 0 
dt 47 r<r(r) dr 2 

for both the hemispheric poloidal flux 4Y(r) ~ 

2nP^ d J d(cos9)B r (r , 9), and the toroidal flux $^(r, 9) ~ 

27ri? w d j{ ivrd dr'B^r' ,9) (defined locally in 9). 

The decay of the flux is shown in Figure [20] along with 
the toroidal magnetic energy. We choose a simple initial 
magnetic configuration with toroidal flux concentrated 
between depth 0.5Ar and Ar, where Ar = i? w d — r(AM) 
is the depth of the magnetized layer. Both the toroidal 
and poloidal fluxes decrease by a factor ~ 1/7 over a 
time interval 10 times longer than the local decay time 
(p/r](Ar), as measured at depth A r. 

At the top of the magnetized layer of both model WDs, 
tdiffjp ~ 10 7 yr. We conclude that the magnetic field, 
although initially buried, will emerge at a moderate age 
in both intermediate-mass and high-mass white dwarfs. 

In the 0.55 M 0 remnant of the solar-mass star tdiff ,l P 
extends up to ~ 1 — 2 x 10 9 yr at the base of the magne¬ 
tized layer; whereas it reaches a maximum ~ 2 x 10 8 yr 
in the 0.87 M 0 remnant of the 5Af 0 star. Although es¬ 
sentially all WD progeny of intermediate-mass stars will 
have strong buried toroidal magnetic fields - due to the 
high angular momentum remaining at the end of the MS 
- some decay in this buried field is expected at ages ex¬ 
ceeding ~ 1 Gyr. 

These ohmic timescales depend on the thickness of the 
C/O layer accumulated on the thermally pulsating AGB. 
The WD masses we have obtained correspond to a nor¬ 
malization ? 7 b = 0.05 to the mass-lo ss rate on the AGB 
using the formula of iBlockeil (|1995l ). There are some 
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Fig. 19. — Similar to Figure flsl but for 0.87Mq WD remnant 
of5 Mq star. Mass processed during the tail of the first thermal 
pulse, which produces the majority of the magnetic flux (Figure [8]) 
is bounded by the dashed green line. 

suggestions of a lower normalization (e.g. r] B ~ 0.01, 
based on the abund ance of Li-rich giants in the LMC: 
lYentura et al.l l2000f) , which implies a larger final C/O 
core. For example, increasing the remnant mass of the 
5 Mq progenitor to ~ 0.90 Mq results in an ohmic dif¬ 
fusion time ~ 4 x 10 8 yr at the base of the magnetized 
layer). One sees from Figure l20l that any toroidal flux 
buried in such a WD will have decreased by a factor 
~ 1/5 from its post-AGB value by an age 3 x 10 9 yr. 

Whether this decay of a buried field corresponds to 
significant decay of the surface dipole magnetic field - or 
even to growth - remains uncertain because the two com¬ 
ponents do not evolve independently. Ohmic diffusion of 
the magnetic held is so slow that the magnetic held easily 
makes a continuing hydromagnetic adjustment to some¬ 
thing close to magnetostatic equilibrium. An interchange 
between the toroidal and poloidal components is possible 
without any change in "H. 

7. SUMMARY AND COMPARISON WITH ALTERNATIVE 
THEORETICAL APPROACHES 

We have investigated how the rotation and magnetism 
of a giant star respond to the inward advection of a small 
amount of angular momentum by deep convective plumes 
in the extended envelope. Some features of the dynamo 
process operating near the core-envelope boundary are 
unique to giants: namely those driven by the intense ra¬ 
diation flux and the rapid inward drift of material to the 
burning shell(s). The buoyancy that is induced by radia¬ 
tive heating of magnetized material strongly enhances 
the rate of mixing between core and envelope. 

Even a weak helical magnetic held enforces a rota¬ 
tional coupling between core and envelope, and main¬ 
tains nearly solid rotation in the core. The inhow of 
mass into the core is accompanied by an inhow of mag¬ 
netic helicity, which is needed to sustain a stable, large- 
scale magnetic held in the white dwarf remnant. The 
compensating helicity is lost through the surface of the 
star. 



(^wd r )/(^wd AM 



t UpAfliM 


Fig. 20.— Top panel: ohmic diffusion of toroidal magnetic field 
intially concentrated at depths 0.5Ar to Ar <C R w d below the 
surface of a cold white dwarf. Here Ar is the maximum depth of 
the magnetized layer of mass AM. Bottom panel: decrease of the 
toroidal (or radial) magnetic flux threading this layer, along with 
the toroidal magnetic energy, both normalized to the initial values. 
This calculation does not allow for an interchange between toroidal 
and poloidal fluxes resulting from a hydromagnetic instability that 
might accompany ohmic drift. 


This strong core-envelope coupling is easily maintained 
in subgiants and core He burning stars, as we found in 
Paper I. It is also sustained near the tips of the RGB and 
AGB in isolated intermediate-mass stars if angular mo¬ 
mentum is distributed broadly within the slowly rotating 
parts of the envelope. We found that the threshold for 
dynamo action is easily reached when H oc r~ 2 at Co < 1. 

Injection of angular momentum from a planetary or 
stellar companion is needed to maintain a magnetic core¬ 
envelope coupling during the giant expansion of late-type 
stars that lose most of their initial spin to magnetized 
winds (corresponding to Mzams ^ 1.3 Mq). 

Although a strong coupling does not depend on angu¬ 
lar momentum injection in intermediate-mass stars, the 
remnant magnetic field may be significantly enhanced by 
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the absorption of a massive planet or brown dwarf, due 
to the relatively low mass of the outer shell that becomes 
magnetized. 

These results are broadly consistent with an increased 
incidence of strong magnetism in i) massive WDs, and 
ii) accreting WDs in CVs, which experienced a tidal in¬ 
teraction before and during a common-envelope phase. 

The strength of the remnant dipole magnetic field is di¬ 
rectly tied to the magnetic hclicity that has accumulated 
in the hydrogen-depleted core. Helicity growth depends 
on a combination of i) efficient angular momentum trans¬ 
port by poloidal magnetic fields through the outer core; 
and ii) a persistent radial angular velocity gradient in the 
tachocline that is driven by angular inhomogeneities in 
the rotation rate in the surrounding convective envelope. 
The helicity flux is roughly proportional to the rotational 
kinetic energy of the inner envelope. 

The WD dipole fields so obtained can exceed 10 7 G 
in both the 0.55 Mq WD remnant of a solar-nrass star 
that absorbs a Jupiter, and the 0.87 Mq remnant of a 
5 Mq star that has no binary interaction (Figure 0. In 
the former case, the engulfment of companion much more 
massive than Jupiter requires exceptional circumstances, 
and this field can be considered a approximate upper 
limit. In the latter case, the engulfment of much more 
massive companion, even ~ 0.1 -Mq, may be significantly 
more common. This results in an increased remnant field, 
reaching as high as ~ 10 8 G (Figure [TUI). 

These field strengths do not approach the strongest 
measured in WDs (~ 10 9 G), which however are very 
rare. The strongest fields may therefore originate in more 
extre me merger events s uch as the collision of two WDs 
(e.g. iGarcfa-Berro et akll2012i) . 

The engulfment of an Earth-mass planet by an evolved 
solar-mass star, which barely augments the angular mo¬ 
mentum of the star, results in a WD dipole field of ~ 10 6 
G. In this case, we find that most of the magnetic helicity 
is deposited during a brief interval at the start of core He 
burning. The star falls below the threshold for dynamo 
activity near the tips of the RGB and AGB. Even a mild 
enhancement of wind-driven angular momentum loss, or 
a softening of the rotation profile in the envelope, would 
eliminate most of this MG magnetic field. 

The magnetic field generated on the giant branches is 
initially buried in the WD (in cases where AGB core 
and envelope decouple before the final thermal pulse) 
but diffuses ohmically outward over the first ~ 10' yr. 
A strong internal toroidal field begins to decay at an 
age ~ 4 x 10 s yr in a ~ 0.9 Mq WD, as compared with 
~ 2 x 10 9 yr in a 0.55 Mq WD. We find that the net 
decay is by a factor ~ 0.2 at an age ~ 3 x 10 9 yr in a 
0.9 Mq WD. The initial decay of the toroidal field may be 
accompanied by transient growth of the external dipole, 
which is possible at constant internal magnetic helicity. 
On the observatio nal side, t he compilation of magnetic 
WDs in iKeoler et all (l2013l f only extends to an age ~ 
1 Gyr, and the incidence of magnetism in older WDs 
remains uncertain. 

We find that the spin angular momentum of the rem¬ 
nant WD of a late-type star can, in some circumstances, 
anti-correlate with the angular momentum of the star 
during the AGB phase. In stars that retain only a mod¬ 
erate angular momentum at the end of the AGB (less 


than the orbital angular momentum of Jupiter), this de¬ 
coupling sets in before the final helium shell flash. The 
final WD spin then depends most directly on the criti¬ 
cal Coriolis parameter below which the envelope dynamo 
fades away in the inner envelope. In the case of absorp¬ 
tion of a Neptune-mass planet during the post-MS evo¬ 
lution of a 1 Mq star, P wc j lies between 0.1 and 1 d. This 
residual spin rate is a lower bound, as it does not include 
a contribution fr om a final anisotropic contraction of the 
envelope (jSpruitl 1199 8). 

The situation changes if the star is massive enough to 
retain most of its natal angular momentum, or if it ab¬ 
sorbs a Jupiter-mass planet or brown dwarf. Then core 
and envelope should remain coupled beyond the final he¬ 
lium shell flash, and spindown of the core can be directly 
related to a transient reduction in core mass. The an¬ 
gular momentum that is retained by the core during the 
contraction of the envelope depends on whether the core 
experiences a late re-expansion due to a delayed helium 
shell instability. A transfer of mass M env e cor eAf* to 
the envelope, combined with ejection of most of this enve¬ 
lope material, implies strong spindown of the core. The 
net result is that the final WD spin period can exceed a 
year in our 1 Mq model. This final spindown was found 
to be significantly smaller in the 5 Mq model, resulting 
in P w d shorter than a day, because of the smoother tran¬ 
sition to the post-AGB phase. 

7.1. Some Outstanding Issues 

Sensitivity to strength of angular-momentum pumping 
in giant envelope. The Kepler asteroseismological data 
give a calibration of the envelope rotation profile in sub¬ 
giants - as investigated in detail in Paper I - but not in 
stars larger than ~ 10 Rq. Here we have considered the 
consequences of uniform dJ/dM in the slowly rotating 
parts of the envelope. This zone with Co < 1 extends 
deep into the star near the tips of the RGB and AGB, 
unless the star interacts with a relatively massive binary 
companion. A rotation profile oc r ~ 2 in the outer 
envelope is supported by the numerical simulations of 
iBrun fc Palaciosl (120091 ) . 

The orbital angular momentum of a Jupiter can be 
compared with the minimum rotational angular momen¬ 
tum that will sustain a dynamo near the tip of the AGB. 
Starting the planet at a* = 1 AU, the excess is about a 
factor ~ 30(Cobenv,crit/0.3) _1 in our 1 Mq model. The 
shallowest rotation profile which could maintain a min¬ 
imal Coriolis parameter Cobenv ~ 0.3 in an envelope of 
depth P-benv / R* ~ 0.01 is fl(r) oc r -4 / 3 . In this situ¬ 
ation, the amplitude of the helicity flux into the core, 
which is proportional to S2 2 env (Equation (l35l) J. would 
be reduced by a factor ~ 10 -3 . The corresponding re¬ 
duction in Hdipoie, by a factor ~ 0.03, would still allow 
for fields of order MG (see Figure 0, but not magnetic 
fields >10 MG. 

A related consideration is the steepness of the angular 
velocity profile in the inner envelope, where Co > 1, and 
its effect on the upper range that is obtained for the WD 
magnetic field. In this paper we did not simply take the 
profile that was determined in Paper I from asteroseismic 
models of sub-giant and helium burning stars. During 
these relatively compact evolutionary phases, an inner 
rotation profile Sl(r) oc r _1 corresponds to Cobenv ~ 10- 
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30. This profile is consistent with the inward pumping 
of angular momentum by deep convective plumes in an 
adiabatic envelope with gravity g(r) oc r — . A lower 
inner Coriolis parameter would be maintained near the 
tips of the RGB and AGB if the same rotation profile 
were sustained there. 

In fact, the rotation profile is expected to steepen 
to O(r) ~ r~ 3 / 2 in the presence of a similar convec¬ 
tive structure, due to the contraction of the core and 
the steepening of the gravity profile. For this reason, 
we allow the inner index in Equation (jTj to increase to 
a = 3/2 during the expanded phase where magnetic he- 
licity growth is concentrated, and the remnant WD spin 
angular momentum is determined. The net effect is to 
raise the upper envelope of the magnetic field distribu¬ 
tion in massive WDs from ~ 3 MG (a = 1) to 10-30 
MG (a = 3/2). The limiting magnetization and spins of 
lower-mass WDs hardly change over this range of a. 

What happens to any core magnetic field left over from 
the MS phase, when exposed to the turbulent motions in 
a slowly rotating convective envelope (Cobenv <C 1)? 
Consider taking the magnetic flux threading a WD and 
spreading it across the core-envelope boundary of the 
progenitor during its AGB expansion. Even a MG mag¬ 
netic held in the WD becomes dynamically insignificant, 
B r (R benv) ~ 10 _4 (47rpUcon) 1 ^ 2 - Such a weak poloidal 
magnetic held would rapidly diffuse over the boundary 
sphere, causing huid elements that are threaded by ra¬ 
dial held of opposing signs to be mixed together. In this 
situation it is difficult to see how the material how from 
the envelope into the core could maintain a large-scale 
poloidal magnetic held unless it were self-consistently 
maintained by a dynamo process. A numerical experi¬ 
ment studying the diffusion and reconnection of such a 
seed held across the spherical boundary between radia¬ 
tive and convective layers is feasible and would provide 
interesting results. 

Relaxation of helical magnetic field to the ‘maximum- 
dipole ’ configuration. The surface dipole magnetic of the 
WD remnant has been estimated assuming that ~ 
~ 'H 1 '' 2 . Other magnetic conhgurations are possible, 
especially those with relatively stronger toroidal helds, 
corresponding to » $ r = (H/^) 1 / 2 l|Braithwaitel 
12009( 1. Ohmic diffusion in WD stars with strong toroidal 
helds may cause some transient growth of the dipole held, 
as the magnetic held relaxes to a more isotropic configu¬ 
ration. 

Dispersion in magnetic field strength in massive WDs. 
Although a greater fraction of massive WDs exhibit 
strong magnetic helds than do their ~ 0.55 — O.6M 0 
cousins, a majority do not. Intermediate-mass stars re¬ 
tain more angular momentum at the end of the MS, 
thereby facilitating a dynamo. This angular momentum 
can, nonetheless, be enhanced by the absorption of a 
massive planet, brown dwarf, or low-mass star, or a tidal 
interaction with a stellar companion. Other effects caus¬ 
ing dispersion in the visible surface magnetic held include 
exchange between toroidal and poloidal magnetic huxes, 
and the time spent by the progenitor star during the 
thermally pulsating AGB phase. 

Emergence of ‘frozen magnetic field generated during 
the pre-MS phase (or during MS core convection)? The 
outer 0.3 M q of the 0.55 Mq WD is processed through 



star ^He core (^o) 


Fig. 21. — Dependence of AGB stellar radius R+ (black line) and 
base of convective envelope -Rbenv (dotted red line) on residual H 
mass. Endpoint of Mzams — 1 M© model. 


a deep convective envelope during the RGB and AGB 
phases. The proportion of this material which experi¬ 
ences rapidly rotating convection, and develops a strong 
magnetic held, depends on the initial placement of plan¬ 
ets around the star. For example, if the closest Neptune 
or Jupiter orbits at a distance ~ 1 AU, then only the 
outer half of the processed layer is strongly magnetized. 

The ohmic diffusion time from the base of this pro¬ 
cessed layer approaches a Hubble time, independent of 
the architecture of the planetary system. From this we 
conclude that the magnetization of all but the very old¬ 
est ~ 0.55-0.6 Mq WDs is not signihcantly influenced by 
a ‘frozen’ core magnetic held. 

The assembly history is very different for the massive 
WD remnants of intermediate-mass stars. Now the bulk 
of the WD mass is processed during core H and He burn¬ 
ing, which leave unaltered the magnetic helicity stored in 
the stellar core. In the case of an isolated ~ 0.9 Mq WD, 
a core magnetic held could emerge ohmically through the 
outer ~ 0.05 — 0.1 Mq that is processed on the thermally 
pulsating AGB, but only with a delay of ~ 2 — 5 x 10 8 yr 
(depending on the normalization of mass loss during this 
phase). A growth in the incidence of magnetism in older 
WDs could therefore point to a contribution of such a 
‘frozen’ held to the visible surface held. (But not exclu¬ 
sively so, given the competing possiblity of dipole held 
growth mediated by the decay of a toroidal held that is 
buried in a thin outer mass shell.) 

Implications of Binary Magnetic WDs. The incidence 
of strong magnetism (exceeding 10 7 G and implying clas¬ 
sification as a ‘polar’) exceeds ~ 3 0% in short -period CVs 
with ages exceeding ~ 1 Gyr (lGansickell2005D . Our calcu¬ 
lations suggest that these helds might decay signihcantly 
in more massive isolated WDs. Long-lived magnetism in 
old binary WDs with periods < 2 hr (below the ‘period 
gap’) could be maintained by the accretion of > 0.1 Mq 
from the companion, which would have the effect of push¬ 
ing down the magnetized layers to a greater depth and a 
longer ohmic timescale. 
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7.2. Competing effects of a collapsing low-mass 
10 Mq) AGB envelope 

The final contraction of an AGB envelope, beginning 
at a photospheric radius ~ lO 2 i? 0 , leaves behind a hy¬ 
drogen layer of mass M enVi f ~ 10~ 4 M 0 . Mass loss slows 
down significantly as the star contracts inside a radius 

i? P AGB ~ Rq (Figure El}- 

It is possible that the process driving mass loss is 
so anisotropic that the collapsing envelope carries with 
it a significant angular momentum, equal and opposite 
to that lost during the final stages of the super-wind 
(iSpruitl Il998h . An upper bound to the angular mo¬ 
mentum thus added to the post-AGB star can be ob¬ 
tained as follows. The collapsed envelope will experi¬ 
ence a bar instability and spread outward into a disk 
if its angular velocity exceeds ~ 0.3(GM w d/i?p AGB ) 1,/2 
at radius R p agb- Any angular momentum in excess of 
Ji nax ~ 0.3M en v,f(GM wd I?pAGB) 1//2 is transported out¬ 
ward by viscous torques in the disk, and is presumably 
lost in a wind. 

Setting the envelope angular momentum to this limit¬ 
ing value, and assuming that it is shared with the rest 
of the star at radius 7 ? p agb, we obtain a minimum spin 
period for the cold WD remnant, 


-P wd > -y—— = 0.4 

^max 


M 


env.f 


10 -4 Mq 


^-RpAGB^ 


hr. 


(57) 

A comparison with Figure [13] shows that somewhat faster 
WD spins can be imparted by rotating envelope con¬ 
traction than would result from a cutoff of the envelope 
dynamo before the final intense phase of mass loss. How¬ 
ever, the fastest WD spins may depend on a merger be¬ 
tween two WDs. 


7.3. Impulsive Dynamo Amplification in a Merger? 

The formation of magnetic WDs in stellar binaries is 
considered by iTout et ahl (1200811 . It is noted that the 
injection of angular momentum from a binary compan¬ 
ion (either a low-mass MS star or a substellar compan¬ 
ion) into a giant star will trigger strong differential rota¬ 
tion in its core and envelope. This differential rotation 
may spark the impulsive growth of the magnetic field in 
th e core , furt her deta ils of which have been examined by 
Wic kramasinghe et al.1 ( 20141) . 

iNordhaus et all (|201ll ) consider the tidal shredding of 
a planetary companion close to the giant core. They 
argue that this material initially forms a disk in nearly 
Keplerian rotation about the core, in which a magnetic 
field is generated by the magnetorotational instability. 
The most promising application of this process is to iso¬ 
lated magnetic WDs, since it may be difficult for a planet 
and a stellar companion (later the mass donor in a CV 
system) to co-exist on similar orbits. 

The contribution made by such an impulsive dynamo 
to the net magnetic helicity stored in the stellar core 
can be compared with the longer term effect examined 
in this paper. When the companion is tidally disrupted 
(either within the convective envelope, or within the rar¬ 
efied mantle surrounding the burning shells), the mass 
M 2 deposited would be re-distributed back through the 


convective envelope on the thermal timescale 


tth 


GM coie M 2 ~ 60Q / M 2 \ (R tide \ 
-Rtide-bAGB \0.1 Mq ) y Rq ) 


yr. (58) 


As envisaged bv lNordhaus et all •: 2 01 ll ), a centrifugally 
supported disk may spread on a shorter timescale than 
(1551) , in which case the internal heat generated in the disk 
is advected around without being radiated back into the 
surrounding stellar envelope. Stud ies o f advecti v e accre- 
tion flows in other contexts (e.g. [Blandford fc Begelmanl 
119991 1 suggest that a fraction of the disk material would 
be accreted deeper into the gravitational potential of the 
degenerate core, which may be enough to transiently in¬ 
crease its rotation. 

Since the energy released by H burning easily suffices 
to unbind the ashes from the degenerate core, any such 
rapidly accreted material must remain H rich. On the 
longer thermal timescale, this nuclear energy can be ra¬ 
diated away, accompanied by a relaxation of the outer 
core to its equilibrium structure. The mass SM nuc pro¬ 
cessed by H and He burning during that relaxation is 
small compared with the total deposited mass M 2 . Given 
that a mass M acc is accreted close to the burning radius 
-Rburn ~ 0.1 Rq, one has 


SM, 


GM a 


-Ms, 


nuc 73 1K2 acc 

^•nuc-^-burn 


'2x 10 -3 


Rbu 


0.1 Rq 


-1 


M ac , 

M 2 


M 2 . (59) 


The normalization is about ~ 10 times larger if H burn¬ 
ing is extinguished, so that only the He-burning shell 
remains active. 

We see that a relatively small mass is incorporated into 
the C/O core before the star adjusts to a new equilibrium 
structure. Given that the magnetic coupling between 
core and envelope is sustained, most of the added angular 
momentum in this new equilibrium state will be stored 
in the outer envelope. 

The large-scale Maxwell stress B r B^I 47t in the outer 
core may easily increase by a factor ~ 10 4 in the im¬ 
mediate aftermath of the tidal disruption, when the ro¬ 
tation frequency increases by > 10 2 . However, as the 
envelope dredges up the injected angular momentum, we 
expect the Maxwell stress to adiabatically adjust down¬ 
ward, since little mass flows into the core. Such a down¬ 
ward adjustment could be avoided if the magnetic field 
lines became decoupled from the envelope, but that does 
not seen consistent with the rapid buoyant motion of 
such a strong magnetic field (see Equation (IT9l) l. 

From this perspective, the impulsive growth of the 
magnetic field following the tidal disruption of a low-mass 
star or substellar companion may leave little permanent 
mark on the magnetization of the C/O core. A merger 
during the early post-MS expansion could induce some 
direct hydrodynamic mixing with the hydrogen-depleted 
core, but only when the accumulated helium mass was 
still well below the final WD mass. 

We argue that the dominant effect of such a merger is, 
instead, to i) push the envelope above the threshold for a 
gentler and more persistent dynamo process in low-mass 
giant stars, Mzams /$ 1-3 M 0 ; and ii) to augment the 
inward flux of magnetic helicity in stars which remain 

































rapidly rotating at the end of the MS. 

The merger of two WDs (or of a WD with the core 
of an evolved companion) provides the main counter¬ 
example to these conclusions. Now hydrogen-depleted 
material in the more massive WD experiences direct hy¬ 
drodynamic mixing with the tidally shredded remnants 
of the companion. The rate of such events appears to be 


relatively small, and may therefore only accomodate the 
most strongly magnetized WDs. 

We would like to thank David Arnett and Peter Gol- 
dreich for conversations. This work was supported by 
NSERC. 


APPENDIX 


TWISTING OF A RADIAL MAGNETIC FIELD 


Here we construct expressions describing the large-scale flow of magnetic hclicity in the evolving core of a giant star. 
The rotation profile and magnetic field are assumed to be axisymmetric, but departures from reflection symmetry 
about the rotation axis are considered. 

Our treatment is simplified by considering a weak radial magnetic field with angular flux profile 

r e 

$ r (0) = / sin O' r 2 B r {6')d0' (Al) 

J 0 

and corresponding potential 

A*(0,r) = —$ r (0). (A2) 

r sin 9 

In what follows, we therefore neglect any radial change in $ r (0). 

A dynamo operating near the convective boundary will generally supply a poloidal field of a fluctuating sign, but 
the contribution to A • B in one hemisphere is insensitive to this sign. The contributions from the two hemispheres 
generally have opposing signs. 

Consider the case where the rotation profile is reflection-symmetric, so that convective motions enforce a certain 
dVt/dr at polar angles 9 and 7r — 9. The evolution of the twist in the core depends on the connectivity of the poloidal 
magnetic field. If both angles 9 and 7r — 9 are connected by the same (axially-symmetric) magnetic flux surface, then 
there is no evolution of the net twist along this flux surface. On the other hand, if the distribution of magnetic flux is 
different across the two hemispheres, then a given flux surface will experience a net differential winding. 

The rate of change of "H inside a spherical boundary of a fixed radius r is 

fyij r 

— = - dSr ■ [c0B + (A • B)v - (A ■ v)B], (A3) 


where the bounding sphere is placed just inside the convective boun dary. 

The electrostatic potential (f> appearing on the right-hand side of (IA3I) is driven by the mean rotation of the star. 
The corresponding electric field E = — ($7 x r) x B/c is sourced by 

<j>(0,r) = ~n(r)<f> r {0). (A4) 

c 

Differential rotation generates A r (and thence helicity) via 


dA r 

~dt 




(A5) 


The electrostatic potential cancels off the term A^cj) ■ (fl x r) in Equation (IA3I) . and we are left with 


<m 

dt 


dH dR henv 
dt dt 


dSA^B^ = —8v r 


dS A^B^. 


(A6) 


Here we have subtracted off the velocity of the convective boundary from the radial flow speed v r of spherical mass 
shells: 


<5v = v — 


dRhenv . 

-;- 1 

dt 


(A7) 


On the AGB, v r is negative and larger in magnitude than dRhenv/dt, but we include this correction for completeness. 

Because the poloidal flux surfaces are, in this treatment, fixed in position, it is useful to re-express the toroidal field 
in terms of a twist angle 4>b'- 


Bc/>{r,9) 


d(j>B 

dr 


r sin9B r . 


(AS) 


The toroidal component of the induction equation is 








which simplifies to 
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d f dcj)B\ d 2 (j)B d 2 (j>B dfl dv r d(j)B 

dt \ dr ) dtdr Vr dr 2 dr dr dr 


in the case (considered her e) t hat <l> r onl y depends on latitude. 
Substituting expression (IA8I) into (1A6I) gives 


dt 


—2n5v r 


d9$ r 


d$ r d(j) B 

d9 dr 


(A10) 


(All) 


which mirrors the standard expression for a static plasma. 

In general the radial magnetic flux profile differs in shape, as well as sign, between the two hemispheres. A complete 
cancellation in the integral over 9 in Equation (IA11I) remains possible if (j> B is independent of 9. That would be the 
case if fluid stresses were to smooth out latitudinal differential rotation within each shell, so that dtt/dr, the source 
for (f>B , were only a function of radius. We see that a net flux of magnetic helicity into the core depends on latitudinal 
differential rotation - differential rotation across flux surfaces - which is a generic consequence of convection. 
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